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■M'lits  ii  vector  tlieorv  for  the  l>ai'kscallerint:  of  rlrelroinapielie  radar  wa\<> 
riie  liasic  leclini(|iie  etn|>loyi'il  in  llie  Miliilion  rcqiiircfl  siiniilaling  llie  vep-la 
ti  mnliiirn.  TIti.s  tmvliiiiri  po.<.'«'!i.«‘.s  an  riecirieal  {wrrnillirilv  that  is  p-iierale«l 
ndoiii  process  and  is  cliaracicri/.cd  liy  a particular  prolialiility  density  function, 
radar  l•ackscattcr  coefficient  is  olitaincd  in  terms  of  the  statistical  character- 
mi  mcdiiini.  \ comparison  of  the  theory  with  rx|i<‘rimcntal  data  is  (liven. 
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hisitihl  i.s  I'ivi-n  into  the  iiiitiirt-  nf  tli'|M>lariza(i«*ii,  l>ul  (‘xplicil  rl■^ull^  for  Ihr  (l<-|io|jiri/^il  li-rnih  an' 
not  olitaiiK'tl  •It  thU  tinw-  Ix'cauM*  of  Ihr  roinplrxily  anil  ilifrirully  of  Ita-  Mjhiiioii.  Somr  of  IIh' 
roiirlusionii  of  thix  work  arr 

^ I.  A throrv  has  Ix'rn  tlrvrlofx'd  for  roiiiputiiii:  llir  likr  |iolari/.ril  (III!  and  VV)  radar  l»a«'k- 
sralirr  riM'ffiririii.s  from  rrrtaiii  ty|M‘s  of  \rfiriation  hy  ii.sini'  a \rrlor  rriiormali/.aiion  a|i|iroa<'li. 

2.  INo  rij^orous  (|iiaiititativr  rompari.soii  of  lln'ory  willi  rxjirriiiirril  was  possihh';  howrvrr. 
({ualitativr  runiparisons  iiidiratr  rra.sonahk-  a|;rrrnK'nt. 

■f.  Althoui^  no  rxpliril  sohition  wax  obtairxd  for  the  ib-pdarizalioti  rornpofW'ittx.  it  wax 
Iramrd  that  one  rau.sr  of  depolarizalioii  is  Ihr  anisotropy  ass<H-ialrd  with  llir  rorrrlaliim  ftinr- 
lion  of  the  dn'k'rlrir  flurtnations. 
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'I'Ih'  uiitliorit)  f«»r  |H'rf(iriiiin}’  llic  Hi»rk  dcMTilMMl  in  llii>  ri‘|)iirt  is  •■•nilaiin'il  in 
l’roj<‘«  l 4 \ in  (icoilrlic,  (iarlnf.'raplin  . ainMn-op^jjiiiic  S<  irinv- 

( \nalvsis  <»f  Radar  ItackscalUT  from  Tfrrain). 

TIu-  (In  ory  tlosrriln-tl  is  the  r«*sull  •»[  an  in  liuiisf  slnity  and  is  an  a|i|iliration  «»f  llir 
work  »»f  olln-rs  in  ihc  fudd  of  »-l«Tlroniatiin‘lii-  wa\«-  siallcrin^  from  ramioin  inrdia. 
riu’  result  is  a inatlii'inatical  model  f«)r  «'aleidalin^  llie  radar  liaeksi-aller  from  eerlain 
lyjM  S of  vefietalion.  The  teelmi(|ue  us4-d  in  the  sidution  is  the  renormali/.alion  formula- 
tion. The  author  wishes  t»)  thank  rrofe.s.sor  Koj;er  II.  I.aiifi  of  the  (ieorj'e  Washim^ton 
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of  fh<‘  T.S.  Army  Kiiffliieer  'ropoj^'aphie  hahoratories  (I'SAKIl.)  for  a-ssistanee  in 
reviewin*;  and  th‘\elopin«  tiu'  tln’ory  de.seriln’d  h«‘rein.  this  task  was  performed  under 
tin-  .supervision  «»f  Mr.  Iternard  11.  .'seheps.  (ihief,  Teehmdo>:y  Development  Itraneh, 
Mr.  Alphonse  Kiser,  (ihii-f,  (ieowraphie  Information  .'systems  Division,  and  Dr.  Kenneth 
K.  Kothe.  Direetor,  (Jeojiraphie  .''eienees  l.ahoratorv.  I he  work  was  under  the  general 
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HACKSCAriKRINC;  OF  KADAK  WAVKS  BY  VKOKTAIKI)  IKHKAIN 


iMKootcrno.N 

PuqicMe.  riiis  n-pi»rl  prcM'iil.s  a \ri  tur  tliciiry  for  ifw  lta<'kM'ull(Tint'  of  i*li*rlro- 
radar  waves  from  vc^cialioii.  Tlir  liasii-  ii’clitiitiiK*  emploved  is  llial  of  simii- 
latino  tlif  vt‘p‘la(ioii  witli  a random  medinm  and  tlien  rali'idatinu  llie  M-alterecI  i-lerlro- 
ma^iietii-  field  from  (lie  random  medium. 

Kaekf'round.  At  (lie  Xriiiv  Kn^ineer  Topographie  l.alioralories.  work  on  (he 

proldem  of  radar  liaekseatterin^  from  natural  terrain  .surfaces  lias  Im-cii  in  progress  for 
many  years.  However,  tliis  work  differs  from  previous  work  iM-eaus)*  it  attempts  to 
model  baek.scatter  from  realistic,  extended  taruets  (terrain)  rather  than  sim|de  *jeo- 
metric  objects  or  surface.s.  After  (he  problem  was  ilefined  and  a literature  review  was 
completed,  the  first  work  calculated  the  radar  bai  kscatter  from  a bare  jiround  surfai'c 
with  only  lar^e  umlulations.  After  this,  the  complexity  of  the  bare  (iround  problem 
was  increa.si‘d  by  addin»  on  small  undulations.'  .Next,  (he  problem  of  .si-atterint;  from  a 
slightly  rou^h  surface  with  a lossy  lavxT  was  cvmsidered.^  However,  the  mathematical 
models  develofX'd  for  the.s«‘  situations  are  not  applicable  to  radar  seatterini!  from  vejie- 
tation  or  cultural  features.  It  was  decided  that  wave  si-atterinu  from  vi-<ji-|j(ion  pre- 
s4-nted  the  most  formiilable  problem  and  therefore  should  Im'  attacked  first.  Tlii^ 
report  presmits  in  detail  the  analysis  performed  in  obtaining  an  a|ipro\iniate  mathe- 
matical moilel  for  radar  .scattering  from  vegetation.  I'lie  work  included;  (I)  a brief 
survey  of  the  available  litera(ur<‘;  (2)  (he  development  of  a model  for  vep'tatioii: 
(2)  the  derivation  (d'  e(piations  for  the  rail  ir  backseatter  eiH-fficients  for  the  model 
developed;  (4)  an  analysis  of  parameter  variations;  and  (.A)  (he  comparison  of  (he 
developed  theory  with  some  experimental  data.  I'hroughout  the  study,  (he  X band 
(tl  (iH7,-l2  (illz)  freipieney  ran^e  was  emphasi/.ed;  however,  (he  developed  solutions 
are  also  ap()lieable  to  certain  other  freijueiicv  rej^oiis.  such  as  K.  and  Ku  bands. 

There  are  two  direct  applications  for  (he  re.sults  of  (his  work.  The  first  applica- 
tion relates  to  the  military  ^eo^raphic  anaivsis  probli'in.  It  is  ho|MMl  that  bv  relating 
radar  backseatter  to  an  a|)propria(e  vegetation  moilel.  one  will  be  able  to  ipiantitalively 
analyze  and/or  classify  a^rieulliiral  crops,  forests,  marshes,  and  other  veijetation  from 
radar  returns;  also  it  may  be  possible  to  eomputc  vcfielalion  parameters  such  as  height, 
density  and  moisture  from  radar  returns.  Tiie  knovvled^e  of  such  parameters  could  be 

* K.  A.  Ilrvcnor.  "Marksi attrrinu  of  KIrrIromagnrtir  Wa^rs  from  a Siirfarr  TompoM'iJ  of  Two  Ty|N’^of  Siirfarr 
UoUKHnfs>.’*  I .S.  Army  KrigimM-r  Tofiographir  LalM>ralohr.s.  K<»f1  Urlvoir.  V A.Trrbniral  Krpfirl  F. TI.-TU-Tl • I. 
Oriolwr  1^71.  \l)  737  673. 

K \ lirvrnor.  *'|{a<  k.'srdtirnng  <»f  FIrrIroniagiirlir  \\av«*>  from  a SlighiU  Koiigh  Siirfarr  with  a Fos.m  l.aArr,” 
\ S.  Arm\  Fnginrrr  Topographir  LalHiralorio.  Fort  Mrlvoir.  \A.  In  hniral  Kr|Mirl  KTI.- 1 K-7 1- 10.  I)r<rriilN-r 
I97F  \l»  \0I3  Hf»3. 
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ill  IIh-  niililan  uiial\>i>  *»!'  Ifrraiii.  Mmi.  a <|iiaiililali\<-  hi<m!<  I i an  In-  ii'<'.l  (<• 
Ill-Ill  ili-\i-ln|i  railar  rri|iiirriiii  iil'  fur  iiiililarv  •>ru«ra|iliir  a|i|ilirali<iii>.  I In-  -<  i mnl 
a|i|iliratiiiii  rrlalr>  (ii  llir  lll■\<'l^|lllll‘lll  ol  railar  iiiiaur  'iiiiiilaliwii*.  \l  |iri->4-iil,  mim'-I 
railar  iiiiaui-  .-'iiiiiilalioiis  ilr\i‘lii|i  >!rayli>iii-  \alur'  lia-^-il  ii|iiiii  a liiiiilril  i|iialilaliti‘  (miIi- 
ji‘rti\i-)  iiiiilrr.'taiiiliii*'  iil*  (hr  iialiin-  of  railar  n-lurii'.  ll  inii-  hail  a i|uaiililalivi-  llirnry 
to  i‘\|ilaiii  (hr  iialiirt-  of  railar  ri-liirii>  dial  i'i>iii|iari-il  ri'aMinahly  Mi  ll  Mith  •■V|M-riiiii‘iit. 
thru  thr  ‘'raytoiir  \alurs  roiilil  Im-  av-iaiinl  morr  ohjrr(i>rly . VUii,  riiiiilaiiiriital  iiiiil'-r 


s(atiiliii|!  Ilf  railar  sr  attrriii"  from  tr»r(alioii  Mill  |>rohahl\  Irail  to  otlirr  }'aiii>  not  yrt 
roiitriii|ila(ril. 

Ih-forr  |)roc«*e<lin>;  to  (hr  ilr\rlo|inirn(  of  a \r}:r(ation  nioilrl.  a hrirf  Mir\ry  of 
past  roiitriliiitioiis  to  (his  priihlrin  Mill  Im-  >9vrn.  (aiiisiilrration  Mill  first  Im-  ui\rn  to 
hoth  hivh  anil  Iom  frripiriu-y  ri-poiis.  I'hr  Iom  fri-ipirnry  rrpon  mIII  roiisi>(  of  thr 
V UK,  I IIK,  anil  I.  Iiamls;  anil  (hr  hij:h  rn-i|nrnr\  rr^oiis  Mill  ronsist  of  tin-  \.  Kn.  K. 
anil  Ka  hanils.  'I'hr  Iom  fri-i|iii-nry  ri-^on  Mas  stnilinl  for  hoth  ronininniration  anil 
rrmoti-  si-nsiii*:  purpos<-s.  Many  rrsi-arrhi-rs  ha\r  rxaniinril  Ma\r  propagation  in  juiij'lt's. 
\inon>:  thrsi-.  I’oinnls  anil  l.a  (Ironr-*  fonnil  that  in  thr  \IIK  n-aion  (Irss  than  2(MI 
Mil/)  thry  ronlil  rrplarr  (hr  jnnair  hy  a lossy  ilii-Irrtrir  slah  havina  thr  a\rra«r  ilirlrr 
trir  ronstant  of  thr  jniialr.  Ilaan  anil  I’arkrr^  allrinptril  to  nii-asiin-  thr  a\rra«r  ih 
rirrtrir  proprrtirs  of  a forrst,  anil  'raiiiir*  anil  I'ainir  anil  Ifrnrr*’  sIiomi-iI  that  (In 
lati-ral  Mavr  plavi-il  an  important  part  in  point  to  point  roinmnniration  in  thr  jnnalr 
for  this  frripirnry  ran»r.  In  aihlition.  an  r\lrnsi\i-  mrasnrrmrnts  prin^ain  Mjsrarrinl 
out  hv  thr  Atlantir  |{rs<-arrh  (lorporation’  in  thr  jnn»lrs  of  Thailanil.  Thr  rr»iil(s  of 
thosi-  stuilirs  anil  othrrs  an-  snmmuri/ril  in  thr  priM-rnhnas  of  thr  “Workshop  on 
Hailio  Svstrms  in  Forrstril  anil/or  Vi-artatril  Knvironinrnts.’'**  In  tin-  I,  haml  n-amn. 
thr  lossy  ilirlrrlrir  slah  mrlhoil  hrraks  iIomii  iM-ransi-  of  inm-asi-il  si-atti-rin^  rffn  t-. 
In  this  rraion.  Kosi-nhauin  anil  UomIi-s’  moilrlnl  a forrst  Midi  a ranilom  ninlinm  ami 
romputril  thr  harksrattrrina  usina  a sin«lr  srattrrina  (hrory.  In  aihlition.  I»n'®  roiii- 

I).  J.  Pound.-,  and  A.  II.  l.a  (ironr.  ('.nnsiileriiyi  h'otett  I rfrelaliiin  at  an  Imprrfrcl  Ihrloririr  Slah,  Klrrlmat  (’.lUi- 
nernru!  Kc^earrh  Lalforatory,  t'nivrrsit>  of  Trxa>.  AuMiii,  Kr|K)rt  Ma> 

^ 4*.  il.  Ilaiin.  anil  M.  I'arkrr,  Frasihiiity  Slutiy  on  thv  t xv  of  Oprn  ii  trr  Transmuxton  (ji/mh  ilori 

to  Measure  thr  F.lrctricai  f*roprrtirs  of  Vegrtation,  Statifuni  UrMcanli  ^‘|lr<ial  I rrliiiiral 

Krport  13,  AiikusI  l%f). 

Thf'fxiorr  Tainir,  “4)ii  Kadio  \\avr  IVopa^alion  in  Fori’>l  KiivirunmcnlM.** /f.f.F.  Tranxat  ttonx  on  .Inlf’nrMtJt  un</ 
l*Topaffntum,  Vn|.  No.  6,  NovrmlM*r 

I).  DrncT,  and  T.  Tamir.  “Radio  l.owi  of  Latrral  \\avr>  in  Forr>|  Knvironmrnl.M/' Radiff  .ViVnrr*.  No.  1. 

' Jan.-iky  and  Rail(‘\ . Tropical  l*ropafiatum  Hcsearch.  Final  Report.  \ oliiinr  I.  Fiitfinrennc  l>r|iartmrnl  of  \llatilif 
Re.srarrh  (lorporation.  Alexandria.  VA.  \IUi4dt3lH. 

o 

J.  14.  Wail.  R.  II.  t)ll.  and  I'.  Telfer.  (F.dilors).  liorkxhop  on  Hnttio  Syxtrms  in  httrcstnl  aiuilor  I •'urtatni 
Tnvironmrnts^  Teehnieal  Report  No.  \(!()- A(3k| -7  f.  Fel>niar\  l*)7l. 

^ S.  Rosenbaum,  and  I,.  Itcmle.-.  "('.liiller  I4e|uni  from  \eiielaled  \n*as.  ‘ U TT.  Tranxactionx  on  tnr('nrHi5  (ind 
t^rofHifiation.  \ ol.  .M'-22.  No,  2.  Mareh  I‘t7l. 

Fi-.|en  Du.  Kaytrifih  Scattering  from  I. eaves.  Seirntifle  Report  No.  I . Thr  Ohio  Slate  I niversit\  Fleet ro  .S*  irnr* 
Laboratory,  21  .ianiiar\  I'tbO. 


(hr  Itavlri^li  M-allrriii;;  fruni  a forrsl  itf  lra\(>  of  arhilrar)  i>ririila(ioii.  \ii 
rxaiiiiiialion  of  (hr  h(rra(iirr  in  (hr  hiuli  frri|uriir\  rr^ioii  hhow>  (ha(  trr\  li((lr  (hrorr- 
(iral  work  ha>  In'rii  <loiir.  \(  (Im-m*  hiulirr  frri|iirin'ir.».  (hr  rffrr(^  «»f  M-a((rriiit>  iNTtdiir 
\rr\  iiii|)or(aii(.  In  far(,  i(  ran  U-  .•.hown  (ha(  (lir  Miiuk-  M-allrriii):  (hror\  of  Kommi- 
Itanni  ami  Itowlo  !.>  no  h>nui‘r  ailrc|nalr  >inrr  innilipir  M-a(lrrinu  rffrr(^  lirroinr 
iin|ior(an(.  UrrrnlU.  Ilonn  ami  I’rakr"  allrn)|i(rii  (o  imnlrl  a fon-M  in  lhl^  frri|iirnr\ 
rr^on.  'I'hr)  look  niiil(i|»lr  M-adrrini:  in(o  arronni  h)  rrilu«;iiit:  (hr  lhirknr<>>  of  (hr 
M'urlalion  ami  Ihrn  iim'iI  a >inu|r  M-allrrin>!  ihrorv  <mi  (hr  rrdnrril  ir|;!rlalion  la\rr. 
Kailar  harkM-allrr  iiiraMirrnirnl.o  in  (hr  hiuh  frr)|nrnr_\  rrifion  ha\r  iN-rn  inadr  from 
hrid^rs.  (riirk  hoom>.  and  airrrafl.  I’rakr,*^  (oNtd\rar  \rro»|iarr.'*  I lah> and 
oihrrs  havr  madr  nirasnrrmrnl>  of  harkM-allrr  rm-ffirirnl.*  a*  a fiinrlion  of  iiiridrnrr 
an»lr  for  a^indlnral  rrop^,  lrrr>.  mar>hr>.  amt  ol'irr  vr^^-laliiHi.  A rriirw  of  iIh- 

>arions  inra>nrrinrn(  |>ro;n'am>  i>  ronlainrd  in  a n-porl  h\  Kiiit’  ami  An 

rxaminalion  oi  (hr  liirrainrr  >how.s  (hal  harkM-allrr  nira>iirrnirnl>  from  forr»l!<  ha\r 
Im-i-ii  madr  h\  M‘\rral  |M-oplr.  In  many  raM‘^.  howr\i-r.  (hr  aiiciilar  ro\rra»r  wa* 
limited  and  drpolari/.ation  informalion  wa>  mil  aiailahlr.  AImi.  (hr  railar  ralihralion.- 
iiM-d  in  (hr  \arion>  iiiraMirrmrnts  |iro^am>  arc  diffrrrnt  and  inronM>lanl.  mi  that  i-om- 
parison  of  (hr  ri-Mills  aiming  in\rsliralor>  i.-  iliffirnit.  In  addition,  (hr  ilala  wai- 
arrompanird  h\  \rr\  lilllr  ipianlitalivr  ^onnd  Iriilh  informalion.  rhi>  i>  parliriilarl) 
Iriir  with  rrspri  l lo  rirririral  paramr(rr>.  stirh  a>  rompk-\  dirirririr  ron>lant!>  of  (hr 
vi-rrlalion. 

In  (hr  nrxl  few  para;n'iipl<>.  wr  will  model  a fon->l  h)  rmploun;!  a random 
mrdinni.  and  wr  will  disriiss  lr<-hni(ptr>  (o  ralrulatr  (hr  radar  UarkM-altrr  rorffirirnt 
from  lliiii  moilrird  forrsl. 

MfMlrlini;  of  N'rvriaird  Terrain.  In  orili-r  lo  ralriilair  (hr  radar  harkM-alli-r 
i-(M’ffirirnl  from  \r};rtal«-d  terrain,  (hr  rlrrlri<-al  |M-rmi-aliihl\ . |M-rniilli\il) . and  ihr 
romlm-liv ity  of  (hr  \r;irlalrd  \oliinir  mn>l  In-  >prrifird.  The  prrmrahilitx  for  mo^l 
forrsis  i.H  i-onslanl  (hroii^houl  and  ran  Ih-  lakrn  lo  have  (hr  Mime  |H-rm(-ahili(y  a.-  (hr 
frrr  spare  \alnr.  Thr  prrmillivil) . on  (hr  olhi-r  haml.  \arir>  a ;n'eat  deal.  >im-r  (hr 
valnrs  for  vr»rlation  are  ipiiti-  diffrrrnt  from  (hat  lor  free  spai-r.  Thr  manner  in  whirh 
(hr  prrmitlivil)  \ari<-s  is  so  romplii-aird  (hat  il  would  hr  a foriimlalilr  (ask  lo  Ir)  and 
olilain  a di-Irrminislii-  spalial  variation  model.  .'>m'h  a di-lrrmini>lir  i-alridalion  woidd 

^ ^ r.  Donii.  jiiil  W.  rrjkf.  - I hr  l.riirrjli/rd  I oinnirll  N-rlitrrr  l.ni'.*  .MtIkmi  iif  rtihaiir  l-.nviriiniiiriil.*'  I nilrd 
National  (lommittn*.  InlrniatMmal  Snrnliftr  KaHio  I ntoti.  >pnnt:  Mri-ttnu.  HI-15  junr. 
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rf(|iiiri’  llic  fxart  loi-atioii  of  cvcrx  It-ai',  Imic.  ami  Irt-t-  Iriiiik  in  a lurr.'l.  K\i-n 

if  (lii>  Mt'rr  |i(>s>ilil«-,  Ilic  ralriilaiion  of  llii*  !>i'all<‘ri'<l  rli‘*'lr<inia“nrli('  fii'lil  from  >iii  li  a 
(-om|iliia(i'<l  ilificrtrir  variation  vvonlii  U*  an  riiorinon>  nmli  rlakniu.  \nil  llirii.  Ili<‘ 
ri'Mill.-  wonlil  Ih-  valid  oiilv  for  llir  |>artii'iilar  forr>l  wlirrt*  llir  ilirl<-rlri<  variation  wa^ 
mi'aMiriMl.  If  radar  liai  kMaltiT  |imlu  lion>  vvrrr  waiitrd  ov«-r  'oiiir  otiii'r  forr>t.  tin* 
rntiri-  rom|dirali*d  drl<‘rinini>lii‘  priH-rdiirr  would  liavi*  lo  Im-  rr|M-al<-d.  In  order  to 
liypa.v.  the  diffieiilties  of  the  di  terniini>tie  Mihllion.  we  will  eoiiMiler  the  dieleetrie 
variation.'  in  a fore>t  to  Im*  ucncrated  h)  a randoni  proee>.'.  The  'Viiihol  to  Im*  im'il  lor 
the  relative  dieleetrie  eon.'lani  will  Im-  (i).  I'lie  veetor  f iinlieati>  that  the  relative 
dieleetrie  eon.'tani  i^  .'onie  arliilrarv  fnnetion  of  three  ninliiallv  |M  r|M-ndienlar  'patial 
eoordinate.-^.  We  will  now  replai  e the  fore>l  with  a rainloin  niedinin  in  whieh  (r)  i' 
«enerated  hv  a randoni  proee.'.'  Iiaviii*:  the  >anie  mean.  and  >landard  deviation,  e^. 
as  the  ori^iial  forest.  To  Im*  .'|M*eifie.  we  will  let 

e,  (r)  = e,  + e.  ^(r)  (I) 

where  is  a stationarv  random  proee.'s  with  zero  mean  and  nnil  varianee.  A simple 
ealeiilation  shows  that 

<e,  (r)>  = e,  Ci) 

<(e,  {t)-ej^>  = e/  (•») 


where  the  hraekets  are  i's<*<l  to  iiidieate  the  priK-ess  of  takiiia  a statistieal  avera'^i*.  In 
addition,  it  will  he  neeessary  to  a.'.snme  some  eorrelation  fnnetion.  Itfr  - l').  for  A<(.r). 
1'he  definition  of  the  eorrelation  finielion  ean  Ih*  slated  as  follows: 

< /i(l)  > = M (r -r')  <d) 

The  eorrelation  fnnetion  has  the  properly  that  it  approaehes  zero  as  Ir-r’l  approaeln*' 
infinity.  I'lie  eorrelation  distanee  (or  eorrelation  len<'th)  is  the  distanee  at  wlmh 
H (x  -x’)  decrea.s<'s  to  e"'  times  its  maximum  value.  When  the  eorrelation  fnnetion  is 
direetion  dependent,  then  the  eorrelation  fimelion  is  considered  lo  Im*  anisotropic.  If 
the  correlation  function  is  independent  of  liireetion.  that  is  it  de|M*nds  only  on  Ir-r  u 
then  it  is  considered  an  isotropic  eorrelation  function.  In  {feneral.  the  correlation 
lentil  is  a measure  of  the  average  size  of  the  particles  in  the  random  medium. 

The  conductivity  in  vegetation  will  also  vary  from  point  lo  point.  However, 
tin*  variations  in  eondiietivily  should  not  pro«lue«*  major  effects  at  tin*  frecpieneies  of 
interest.  The  major  effect  of  eondiietivily  is  its  average  value,  o^.  which  helps  alli'ii- 
uate  the  waves. 


Ut'Utro  to  itiM'iiNi  M>mt*  ol  lh«‘  lliat  < an  !«•  iimmI  Io  cal 

fiilalf  lln‘  liarkst  allcr  ctM'ltiriciit,  wr  >liall  t-oin|)iil<‘  .sonic  l\|iical  values  of  o^,  anil 
lor  a ilcciiliioiis  torcsl,  to  ilclcrniinc  llic  orilcr  of  nia^nilmlc  of  IIicm-  |iaranictcrs. 
We  will  a.ssiiiiie  a loresi  made  ii|i  of  leaves  alone.  Kaeli  leaf  will  he  eoiisidereil  Io  he  a 
reelaii^iilar  solid  in  sliape.  with  the  dimensions  of  w .\  w \ I.  'I'he  dimensions  of  the 
entire  lores!  will  he  I,  \ L \ I).  The  spaeinu  iM'tween  leaves  in  the  hori/ontal  direetion 
is  d;  while  the  spaeinj'  in  the  vertical  direelion  is  a.  The  .'elalive  mean  eom|ile\  dielee- 
trie,  e^,  can  then  he  eom|iuled  as 


where 


\ 

\ 


V 


A 

fc 


c 


A 

e = 


' v \ \ f A 


= e - j (o  /caJ  e ) 

d J ' d o' 


volume  of  one  leaf  = w \ vv  x I. 

total  luimher  of  leaves  in  the  forest. 

relative  complex  dielectric  constant  of  one  leaf. 


(•■») 


V ^ = volume  of  the  forest  which  is  solely  air. 

e ^ - relative  dieleelrie  eonstanl  of  air  = 1 .1) 

\ = the  total  volume  of  the  forest  = I.  x h x I). 

u>  - radian  frerpiem-v. 

= permittivity  of  free  spar-e. 

I'imlinj'  the  <lieleelrie  constant  of  leaves  of  any  type  is  diffienlt  to  ilo.  since  .so  little 

data  has  her  n ohtained.  One  »»f  the  f«-w  reliahle  sources  availal  '•  vv.as  used.  i.r-.  Itroad- 

hurst'*  who  im-asured  the  complex  permittivity  of  tulip  tree  leaves  and  hram-hes  for 
a hroad  ran^e  of  freipiencies.  I’rom  his  data,  w«:  find  that  in  the  .\-hand  re<rion.  the 
relative  complex  permittivity  of  tulip  tree  leaves  is  approximately  TO-jlO.  I'sinp  this 

valur;  alonj;  with  tin’  foMowin*;  parameters,  we  can  compute  t . 

* « 

Let 


L = 1.“)  meUrrs 

I)  = d imders 

w = 4 '•entimeir’rs 

t = (t.l)l  centimeter 

d = 4 c(!ntimelers 

a = ()..'}  meter 
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'riu‘  aln»v«*  sc‘l  ol'  iiuiiiIkts  will  allow  l'<»r  llic  calnilalion  of  \ ^ \ ^ , aiol  \ j . W lirii 
tli(‘S4‘  paraiiicU-rs  an*  (-alciilatcd  and  itlaccd  in  tli<‘  ahovf  i‘)|nation!<.  wi-  ohiain  for  i;^ 


= I.OOd  - j (MIOOH 

In  a similar  inannrr,  tlir  \ariaiu'i',  6^^  . of  tlir  rral  portion  of  llir  coniplcx  |M'rmilti\it> 
of  tlu‘  forest  ofleavt'sean  l)e  eompnted  l)y  nsin^  tlie  following  eipiation: 


e'  = 


V\,  (e/  {e ^ -c^)^ 


{(>) 


If  the  tiata  ^iv<'n  above  is  used,  can  he  easily  ealenlated  to  ^i\e  the  tollowin^  value: 

e?  = O.ll’f) 


'I'he  ealenlated  value  <>f  aj^'ees  with  tin'  opeti  wir«‘  transmission  line  results  pres<‘nt<'d 
hy  llaj^i  anil  I’arker.'^  It  ean  he  M-en  that  even  thou;'h  e^,’  is  very  lar>;e.  the  a\era«e 
dieleetrie  is  ap|iro\imateiy  one  heeausi'  of  the  small  pereentajie  of  the  total  volume 
that  the  leaves  oeeupy. 


A forested  terrain  will  Im-  modeled  hy  a flat  earth  eovered  hy  a slab  of 
random  medium.  The  random  medium  is  deserihed  hy  its  dieleetrie  variation  fjiven  hy 
equation  ( I ).  A plane  wave  from  free  spaee  is  a.s.sumed  to  he  incident  upon  the  slab  at 
anijle  0.^  with  respect  to  the  normal.  The  objective  tlieii  is  to  calculate  the  iiackscatter 
coefficienLs  for  huri/ontal  and  vertii'al  polarizations.  This  model  of  the  forest  can  Ire 
simplified  further  when  the  frequencies  of  interest  an-  \-hand  and  above.  Ifeeause  of 
the  averafte  loss  due  to  water  eonteni  and  los.s«'s  due  to  multipli'  seatt«'rin»,  the  inci- 
dent wave  will  Im-  completely  ahsorlred  in  most  easi's  before  reaching  the  {iround.  As  a 
result,  we  can  nej'lect  the  effect  of  the  ^ound  and  assume  that  the  backseat  ter  comes 
.solely  from  the  forest  medium.  This  is  <'(|uivalent  then  to  solvin«  the  problem  of  a 
plane  wave  incident  upon  a half  space  of  lo.s.sy  random  media.  .A  nuinher  of  teehniqiies 
will  1m;  briefly  diseu.s.sed  that  could  he  u.s»'d  to  calculate  the  hackseatter  from  a lossy 
random  half  space. 

Single  Scattering  — Horn  Approximation.  This  technique  was  one  of  the 
first  to  he  u.sed  for  the  calcidation  of  scattering  from  random  media,  basically , the 
method  entails  dividing  the  scattering  medium  into  a large  number  of  individual 
scattr-ring  el(;ments.  'I’he  .scattering  from  each  element  is  computed  separately  and 
independently  from  all  other  elements.  Then,  the  total  scattered  field  is  obtained  hy 
a<lding  the  fields  due  to  each  individual  .scattering  element.  Mathematically,  the 

(J.  II.  lUgn  and  II.  W . Parker,  h’easibility  Study  on  Ihe  ll»e  of  Open-Hire  Trantmistion  l.inet.  Capacitor!  and 
(^viiies  to  Measure  the  Electricat  Properties  of  I'egetatinn,  Stanford  Rrsoarrh  Institute,  S|M*fiaI  Technical 
Re[M)rt  13  August  l%6. 
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pr<>l)l«'iii  involves  an  iiitc^ul  *'(|iiutioii  I'or  llie  scaltereii  field,  in  terms  of  tlie  total  field 
(ineident  field  phis  seattiTed  fiidd).  Tli<‘  seattered  field  »illiin  the  inte^al  is  set  eipjal 
to  zero  in  ordiT  to  oltlain  a first  a|ipro\imation.  'I'lien  the  ineident  wave  aets  like  a 
souree  for  the  seattered  wave.  When  the  seatterinj:  inedinin  is  lossy,  Ihosi*  elements  far 
enon^i  from  the  snrfaee  s«‘e  an  ineident  field  that  is  highly  atlemialed.  As  a result, 
their  eontrihntion  to  the  haek.seattered  field  is  small,  for  this  reason,  the  \evetation 
layer  app«‘ars  to  have  an  equivalent  thiekness  from  whieh  most  of  the  seatterin;;  oeeurs. 

.'since  this  teehnique  does  not  eon.sider  multiple  seallerin;'  effeets  and  is  limited  to 
media  that  are  slightly  random,  that  is  I,  it  was  decided  not  to  explore  this  tech- 
niqne  any  further. 

Multiple  Scattering.  In  the  single  seatterin*'  method  just  di.seiis.sed.  the 
random  medium  was  broken  up  into  many  mutually  independent  .seatterers.  The 
a.ssuinption  of  mutual  independence  amon»  seatterers  is  in  general  not  valid,  and  one 
must  consider  the  coupling  hetween  .seatterers.  This  mutual  coupling  i-an  he  accounted 
for  hy  including  in  the  haekseattered  field  eontrihutions  from  waves  that  have  Im'I'ii 
.scattered  more  than  once.  'I’liere  presimlly  exists  many  methods  that  can  he  usiul  for 
ealeulatinv  the  effi'et  of  multiple  scatterin*».  The  accuracy  of  each  method  and  its 
rei^on  of  validity  are  ipiestions  of  current  resi’arch  and  have  not  heen  completely 
resolved  to  date.  A brief  outline  of  each  method  will  he  pven  in  the  followiii}:  para- 
jiraphs. 

Scalar  Renonnalization.  riie  total  field  in  the  random  medium  is  s«'t  ecpial 
to  the  .sum  of  a coherent  mean  wave  and  a seattered  wave.  'I'wo  (-(piations  are  derived 
from  the  scalar  wave  equation  that  allow  for  approximate  solutions  of  the  mean  wav<‘ 
and  the  scattered  wave.  Thesr*  two  e(|uatiuns  an^  referred  to  as  the  renormalized  ecpia- 
tions  for  the  mean  and  seattered  waves.  'Flu’  .solution  for  the  .scattered  wave  involva’s 
the  mean  wave,  and  so  the  mean  wave  must  he  calculated  first.  In  fact,  the  mean  wave 

acts  like  a source  term  that  "cnerates  the  .scattered  wave.  'Fhe  .solution  of  the  eipiations  . 

for  the  mean  wave  will  brad  to  a dispersion  (‘(piation.  'Phis  dis[)ersion  ecpiation.  when  j 

solved,  will  yield  an  effective  dielectric  constant  for  the  mean  wave.  This  effective  ^ 

dielectric  constant  will  he  complex  eveti  if  the  mean  conductivity  in  the  random 
imulium  is  zero.  The  interpretation  that  has  Ix^en  placaul  oti  this  is  that  the  ima;;inary  I 

part  of  the  effective  dielectric  constant  arises  from  multiple  seatterin*;  effects.  Thus.  j 

multiple  scattering  is  eotisidenul  when  the  .scattered  wave  is  calculati'd.  The  unknown  ; 

factor  in  this  technique  is  how  much  multiple  .scatlcriti*;  is  really  considered.  .Also, 
the  scalar  approach  cannot  allow  polarization  cffci-ts  to  he  considered.  In  order  to 

considt^r  polarization  effects,  one  would  have  to  ii.se  a vector  wave  (‘(juation  and  he  f 

certain  that  the  real  causes  of  depolarization  an-  included  in  the  solution.  ^ 


Vector  Kenomialization.  Fhe  approach  taken  here  is  similar  to  the  scalar 
renormalization  approach,  except  that  now  a vector  wave  ci|uation  is  considcrc<l  atid 


ili‘p<>luri/.ati(>n  terms  run  Im-  ul>luiiie<l.  \s  with  tlie  seular  approaeli.  a (lis|M‘rsi<*ii  eipia 
tioii  is  ohtaim-d,  and  an  eft'eetiv<'  dieleetrie  eonslant  is  Imind.  .Num,  lioMex-r.  the 
effective  dieleetrie  eonstani  turns  out  to  lie  a d liv  d matrix  (tensor):  whereas  in  the 
j scalar  ease,  it  is  a scalar.  Ilasieally,  this  |ii)ints  out  that  we  are  replacing  the  random 

medium  hy  a deterministic  anisotropic  medinm.  The  eletiH'iits  <if  th*-  effective  dielee- 
I trie  tensor  are  found  to  involve  the  correlation  fnnelion  of  the  random  medinm.  The 

I mean  wave  is  then  computed,  as  in  the  scalar  case,  hv  tisinu  the  solution  fri>m  the 

I dispersion  etpiation.  'I'lie  .scattered  field  can  then  Ih-  eali  ulated  hv  nsin^  the  solution 

for  the  mean  wave.  This  techniipie  was  the  one  chos<‘n  for  thoronjih  anaivsis.  ami  it  is 
I the  main  snhjeet  <d'  the  rest  of  this  report. 

Other  Techniques.  A nundicr  of  other  tei  hnicpies  «!xist  that  could  Im-  applied 
to  the  vefjetation  prohlem.  One  such  techni(|ue  is  the  diffu.sion  imdhod.  I’liis  ttiethod 
involves  derivin<i  an  eipiation  for  the  prohahility  d«-n.sity  function  of  the  field’s  ampli- 
^ tude  and  phas<‘.  Once  an  approximate  solution  to  this  ccpiation  has  been  ohtanieil.  the 

ratlar  hackscattcr  coefficients  can  he  readily  computed.  This  techni(|ue  has  ln-eti 
successfully  carried  out  in  the  one  <lim('n.sional  ca.se,'*  i.e.,  the  case  where  it  is  as.sume,d 
that  the  dielectric  constant  varies  only  in  one  direction. 

^ Two  other  techniques  that  could  he  u.s«-d  an'  the  radiative  transfer  method 

and  the  discrete  scatterer  techni(|ue.  The  radiative  transfer  techni(|ue  is  a phenomeno- 
lo^cal  theory  that  attempts  to  calculate  intensity  usinjj  the  a.ssuiiiption  that  the  power 
radiated  from  individual  scattering  regions  adds  incoherently.  Alternatively,  the 
discrete  scatterer  technique  makes  use  of  the  far  fiedd  radiation  characteristics  of  indi- 
[ vidual  scatterers  to  obtain  hackscattcr  characteristics  for  a collection  of  thc.s«'  .seatterers. 

Another  technique  is  computer  simulation  in  which  the  .scatti'red  power 
would  he  computed  hy  numerical  methods  for  many  diffenmt  sample  values  of  the 
random  media.  The  mean  value  of  the  scattered  pow«'r  for  the  ensemble  could  then 
1m'  computed  from  the  individual  sample  values. 

ANALYSIS 

Tlu!  vector  renormalization  approach  cho.sc'ii  for  analysis  in  this  report  recpiircs 
i defining  the  geometry  of  the  scattering  prohlem  and  the  approximati'  solution  of  wave 

propagation  in  a random  medium  of  infinite  extent  in  order  to  derive  ('(piations  for  the 
calcidalion  of  horizontal  and  vertical  polarization  haekseatter  coefficients,  i'.quations 
( 1 2'))  and  ( I 75)  are  the  final  residls  of  all  the  derivations  in  the  following  paragraphs, 
(ionsider  the  problem  of  a plane  electromagnetic  wave  ii\eidcnt  ohliipiely  at  an  angle 
0.  from  a free  space  medium  (z>  0)  onto  a lossy  random  medium  (z  < 0).  I’hc  geonn- 

I H 

R,  H.  Lan^,  “Probability  Density  Function  and  Moments  of  the  Field  Inside  a Dne-DirnenMonal  Random 
Medium,”  Journal  of  Mathematical  l*hysics^  Vol.  II.  No.  1 2,  December  197H, 


H 


e(jl)  = Cq  c^Il)  = £o  t^a  ''' 


Figure  1.  Scattering  Geometry 

Tin;  lossy  random  medium  is  eliaracteriztul  by  a dielectric,  e (r),  wbieli  is  composed  of 
lli(‘  sum  of  a random  part  and  a nonrandom  part.  The  nonrandom  part,  e e , is  equal 
to  the  mean  of  e (r).  The  random  portion  comes  from  p (r),  which  is  generated  hy  a 
stationary  random  process  with  zero  im  an  and  a variance  equal  to  one.  The  objective 
is  now  to  obtain  an  approximate  expression  for  the  backscattered  far  field  in  the  upper 
medium.  This  field  can  then  be  used  to  compute  the  radar  backseatter  coefficients. 
In  ordt'r  to  conqjute  the  scattiTcd  fielil  in  the  upper  medium,  one  has  to  first  obtain 
an  understanding  of  the  nature  of  wave  propagation  and  scattering  in  a random 
tnediiim  of  infinite  extent.  This  will  be  done  using  the  renormalization  formulation. 
The  method  of  renormalization  was  developed  in  detail  by  'I'atarskii  and  ('icrtsen.shteiii.'’ 

The  elecirie  (K)  and  magnetii-  (M)  fields  in  a random  medium  can  be  rclalt'd  by 
Maxwell's  ecpjalions  in  (iim;  harmonic-  form.  The  time  harmonic  variation  will  be 

I.  Talarskii.  anil  \l.  K.  Gi-rlsi'iishli-in.  "FropaKationof  Waves  ina  Mriliiini  with  .Strong  Klurluation.>iof  the 
Krfrarlivf  Index,'*  SonX  I*hysics  J V'ol.  17,  No.  2,  Aunust 
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exp  will  asMinic  a riclaii'iiilar  Carloiaii  coonliiiaU-  syslcm  lo  Ik-  iiiilM'dilcd  in 

III*'  raiidoiii  iiicdiiun  of  iiii'iiiilc  rxlt'tit. 


V X i;  = 

(7) 

V X ii  = 0^  li+  ju;e(r)ii 

(H) 

= c., 

Till'  lino  iiiidiT  a syiidtol  will  lif  usi-d  to  n’|)rrs«Mil  a vfclor  (|iiati(it\. 
of  (7)  and  using  (H),  oiu'  has  llic  following  result: 

'Faking  the  curl 

V xVx  K-k^  K = k„^  c^M(r)  K 

where 

k ^ = CO^/i  € 

F,({uation  (9)  ean  also  he  written  in  operator  notation  as  follows: 

II 

(10) 

III 

<1 

<1 

X 

1 

7“ 

We  will  now  assume  that  the  total  eleetrie  field  (K)  in  the  random 
written  as  the  sum  of  a mean  wave  and  a scattered  wave. 

medium  ean  Im' 

K = < K>+ F; 

(II) 

< K > is  the  mean  wave 

K is  the  scattered  wave 

Obviously  the  mean  of  is  zero.  However,  the  tnean  of  • K_,  will  not  Iw  zero.  1 he 
star  indicates  taking  the  eom|)lex  eotijugate.  Placing  e(|ualion  ( 1 1 ) into  ( 1 0)  and  taking 
the  average  of  the  resultant  expri'.ssion  will  give 

£<F>  = <£h',  > sin<e<£>=  II 

(12) 

10 


\\  lifii  ( 1 2)  i.s  s(il)lra('t(’<l  from  ( 10)  wo  lia\«- 

f li;-<  K,  > I = ^ K-<^  l^,> 

JCK  K>  * ^ K -<£  !•> 

E,  = i:-'  1J<E>  1 + UK,  - <^]L>\ 


(IH) 

(U) 

(ir.) 


riu-  iiivcrsf  ojuTalor  £'*  can  In-  written  in  terms  of  an  infinite  spare  il\ailie  (ireen's 
funetion. 


£-'  ( ) = / ( ) • r<r.r')  .Ir' 

V ' 

where  ilr'  = iIx'.In'iI/.' 


(If)) 


X ) if^  the  infinite  spaee  dyailie  (ireen’s  funetion,  where  the  ilonhie  line  miller^  is 
use.l  to  in.lieate  a dyadie.  'I’he  form  of  T (l,  r')  is  as  follows: 


(IT) 


w-here  J_is  the  unit  dyad.  . A i;ood  ilerivation  of  the  dvadie  (in'.m's  funetion  ean  Im- 
found  in  the  monoiiraph  hy  (ihen-'l'o  Tai.“  The  intejiration  in  (16)  is  to  U-  earried 
out  over  the  entire  volume  of  the  random  medium,  whieh  in  this  ea.se  is  over  all  spaee. 
When  the  equation  for  K,  as  .^ven  hy  (15)  is  plaeed  in  e(|uation  (12).  the  f.)llowin>i 
result  is  ol.tained: 


lX-<f£-U>|  <Ji>=<^X-'  UK,-<fK,>|> 


(IH) 


For  a first  ap[)roxirnation  of  tlie  mean  wavi;,  we  will  eonsider  the  ri<ihl  hand  side  of 
( IH)  to  l)e  zero. 


1£-<JX'U>I  < K>=  0 


(I'n 


An  e.piation  lor  the  seattered  wave  in  terms  of  the  mi'an  wave  ean  he  obtained  l.y 
usinji  (1.5).  The  s«'eond  term  on  the  ri^ht  hand  side  of  (15)  will  he  eonsidered  small  in 
eomparison  witli  (he  first  term  involving  (he  meat)  wave. 


K =X-'  U<K>| 


(20) 


20 


('.hi-n-To-Tai,  Dyadic  (ircpn\  f'linctioiu  in  ElectTomaiinctic  Theory,  International  Textbook  ('.ompany,  1071. 


II 


I (l*>)  ami  (Jll)  an-  llic  |iair  ol'  rrmiritiali/.aliDii  ('(|ualiuii>.  Tlir  i»li><T\a 

lion  llial  "lioiilil  U'  iiia)l<'  i>  lliat  llir  iiiraii  \\a\i'  iIim-s  iiol  |iro|ia»at<‘  with  tlir  |iro|)a;ia 
lion  t'oii.slaiil  k.  Uni-  nni't  I'lr'l  ilrrivr  a ili'|i<'r>ion  i'i|nalion  from  ( I*))  ami  tlicn  .-oUc 
Ihi-  (-ijnalion  inr  jii  rfl'ri  hM-  propagation  i iin-liint.  ’I'hi>  rl'fcclivr  propagation  ron- 
-lanl  will  ili-M  rihi-  lln-  alli  inialion  ami  plniM-  rliararlt'ri>tir>  of  Ilit-  mi-an  wa\r.  Tin- 
M'l'oml  imporlani  oliM-rvalion  llial  >lionl<l  In-  niaih-  i>  that  llir  mran  wave  ai't.s  liki*  a 
Miiirn-  lor  llir  srattrml  wa\r  a^  >liown  h\  ri|nalion  (20).  Wlirn  llir  (|uantilir.>i  for 
f,  ami  £■ ' arr  plarril  in  ( 10),  llir  following  r\prr.s.'>ion  rr-nlls: 

Vx  <K(r)>  - k^<K(r)>  - k^;*  y<p(r)p(r')><K(r')>  • r(r.x')  lir'  =0  (21) 

V ’ 

llir  (|iiantity  '<^/i(r)p(r’)>  wa.i  ilrlinnl  rarlirr  a.-  llir  rorrriation  fiinrlioii.  .M.'O.  il 
.-'lioiilil  hr  notirrd  that  within  ihr  volntiir  inlr<;ral  ihr  mran  vvavr  is  a fnprtion  of  r': 
whrrras  (hr  two  Irrms  iiivoKirii;  Ihr  iiiraii  wa\r  onisidr  Ihr  iiitr<n’al  arr  fiiiirlioiis  of i. 
This  is  iM-raiisi-  rvrrylhina  to  (hr  rialil  of  Ihr  £'•  oprralor  nnist  Ik-  inrlndrd  in  (hr 
intrijral.  ,\lso,  it  ran  hr  shown  that  hrraii.sr  of  Ihr  symmrtriral  |iro|M‘rtirs  of  Ihr 
dyadir  (irrrn's  fiiiirtion  Ihr  following  rqnation  ran  hr  vsrittrn: 

<K(r')>  •r(r.L)=r(r.  r')  •<E(r')>  (22) 

\\r  will  srrk  plant'  wave  solutions  It)  (21)  whirli  will  have  (hr  followin*;  form: 

<K(r)>  = Ar-t!^-t  (23) 

A i*  ronstaiil  vector 
K is  the  effective  (iropawation  ronstaiil 

The  effective  propagation  constant  Jv  is  complex  and  will  1k'  complex  even  if  k is  real 
(o^  = 0).  The  complex  portion  of  K will  attenuate  (hr  mean  wave.  This  attenuation 
has  heen  interpreted  as  arising  only  from  multiple  scaltrriiig  effects,  when  is  zero. 
In  general  then,  there  are  two  factors  that  contrihute  to  the  attenuation  of  the  mean 
wave.  'I'he  first  eontriliution  is  due  to  the  average  lo.ss  in  the  random  medium  that 
ari.se.s  from  a finite  value  of  a.  . 'I'lie  second  eontriliution  comes  from  the  effects  of 
multifile  scattering.  The  mean  wave  that  appears  inside  the  volume  integral  can  be 
written  as 


<E(r')>  = A e-jt'-t 

(24) 

< K(r')>  = A 

(25) 

<K(i:')>  = <K(r)> 

(2fi) 

)2 

J 


\\  ilh  Uu'  aid  of  (22)  anti  (26),  fiination  (2 1 ) can  lx*  written  a.s  follows: 


Vx  Vx  <K(r)>  - k2<EC0>  - K j '*(1-^')  nx-i')'  <t(r)>  ' ilr'=  o 

V * 

[VxVx  ej*  j <ir']  <F.(r)>  = 0 


(27) 

(2H) 


Now',  let  us  consider  an  eleetromapietie  wave  propaxahnu  in  a deterministic 
anisotropic  medium.  'I'lie  electric  field  for  tins  wave  will  be  desitrnuted  by  b„-  I be 
relative  complex  dielectric  tensor  will  be  £.  'I'lie  vector  wave  etpialion  for  tbe  electric 
field  propaj^ating  in  the  anisotropic  metlinm  is  as  follows: 

(VxV  x-k,2  g)K„  = 0 (29) 

where  in  f'eneral  e.can  be  written  as  a 3x3  matrix. 


^xx 

"x> 

Sx 

Sv 

^yx 

.^xx 

^xx 

When  etpiation  (211)  is  compareil  with  cipiation  (29),  we  see  that  they  arc  »t  the  same 
form  and  can  be  made  etpial  if  we  allow  the  following  two  relations  to  hold: 

K - <K(r)> 


and 

i = i + f I{(r-r')  r(r-r')  e iK-<r-i  )dr' 

_ k ^ ~ O S J 

O v' 

riins.  it  can  be  seen  that  the  mean  wave  in  a random  medium  can  be  considered  as  a 
wave  prupa^ating  in  a deterministic  anisotro|)ic  medium.  'I’he  relative  complex  dich'c- 
Iric  tensor  can  be  .seen  to  be  a function  of  the  correlation  function.  |{(r  - r').  <>f  tlii' 
random  medium.  'I’he  dyadic  rireen's  function  can  easily  lx;  written  in  reetanpdar 
eom[)onent  form  so  that  the  individual  components  of  f can  be  determined. 
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r(r-r')=la^i,,  +a^a^  -Ha,aJf,(K)+  |^a,a,  + a^a^  (A-1').(> 7')  + 


^.^x  — 

' ' + a a / + a a 

V\  here 

^x 

is  a unit  vector  in  the  x direction 

is  a unit  vector  in  the  y <lirection 

a-. 

is  a unit  vector  in  the  z direction 

K = 

Ir-r'l  = V(x-x')^  +(y-y')^  + (/. 

■'-'■y  1 f 


(H) 


f2(R)  = (kMi^  -jkU-l)e-J‘^«/(47rkM{^) 


f3(R)  - (3jkl{ +3-k2K^)c-j'‘«/(47rk2K^) 

llie  individual  clt-mcnts  of  the  relative  eomplex  dielectric  teiii^or,  |.  can  now  Ik-  written 
and  are  ujven  Iwdow  : 

^xx  = •'"/K,'  +1^0'  e//  'Mr-r')  |f2(H)+  f3(ll)|  ej!i-(^-i'>dr' 

v' 

= K"  j j f3(H)ej!^-<'-^>  dr' 

v'  ^ ' 

^x/.  = J 'Hf-x')  j dr' 

V ' 

^X  IHr-r')  f3(K)ej'^--(t-t)  dr' 

V ' ^ 

^yy  = '"'/''o'  + l"o'  j >'2(K)+  f3(H)  j ej'^-<K)  dr' 

V.  " “"o'  ^'/  '5(w')  f3(R)«J*^-"<^-"'>  dr' 
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= K'  f//  'Hr-r')  j 

V ' 

(/.-/)  (x-x')( 
1(2  ( 

f3(K),J!i -d-T)  .b' 

= f “(i-i') 

1 (/--•/-')  (y-y')  ( 

j ( 

f,(K)ej*^-<i-^>.|r' 

= kVk^,^  ^k,,^  e/  J l$(r-r')  j fjfH) 

+ f3(R)j 

V 


it  ran  fa.sily  l)*‘  from  the  ul*ov<‘  tliat  lli<*  (licicciric  triisor  e Ik  .syininct- 

ric.  Kai’li  t»l'  tli<*  elcnieiils  of  tin-  dielectric  tensor  can  Iw  s»*eii  to  !>«•  a funeti<»n  of  K. 
'riierei'ore,  not  oiiiy  tiie  iiia^iitiide  but  also  tlie  direction  of  tin*  effective  propat'ation 
constant  is  important  in  determinin;'  each  of  the  elements  in  the  dielectric  t<‘iisor.  W'e 
will  now  attempt  to  obtain  a dispersion  <‘(|uation  from  (20),  which  will  allow  a solu- 
tion for  K-  The  operator  V ^ !’<•  ''rilh*n  in  matrix  form  as  };iv«‘n  below ; 


92 

92 

d/}  dy^/ 

9y9x 

9x9/ 

+ AM 

92 

9x9y 

\9/^  9x^/ 

9/9y 

9^ 

9^ 

/ 9^ 

9x9/, 

9v9/ 

Ux^ 

Hv  iisin*!;  the  above  definition  and  lettiii"  K - E (l)  ^ - A exp  (-jK  • r),  e(piation 
(20)  beeoiTU’s  ° 


K.^  + k 2 -k  2 e - K k -k,^  e. 


z y o XX 


k k - k ^ e 

XV  o V X 


KK.-K"  ^zx 


XV  O XV 


k 2 


2 

o 'vv 


- k k - k ^ e 

V /.  o /.> 


k k - k 2 e 

X /.  O XI 


k k - k 2 e 

y ^ o \ 7. 


k 2 k ^ -k  2 e 

y X O 7.7. 


= 0 (30) 


where  kj^,  k^, , and  k,^  are  the  components  of  k alonf;  the  x,  y,  and  •/.  axes  respective! v. 
bikewist;  A^,  and  A,^  are  the  components  of  A.  It  <'an  he  s*‘eii  tliat  the  system  of 
partial  differential  eipiations  has  been  replaced  by  a system  of  liiicar  alfjebraie  eipia- 
tic'iis.  We  will  now  switch  from  reelan<ruli>r  to  spherical  coordinates  in  our  docription 
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1)1'  dll'  K vi'i'tor.  (In'  ili-l'iiiitioii.'  of  0 ami  0.  tin*  >|)lii-riral  aiialr.'  that 

will  1)1*  ii.«<*(l  to  llrsiTiln*  K. 


z 


'I’lie  an};le  6 is  the  aiif'le  between  the  z axis  and  k.  The  annie  0 is  the  annie  between  the 
X axis  and  the  projection  of  K onto  the  xy  plane.  The  reetaiifadar  eoinponents  of  K 
can  now  be  written  as  follows: 

= K sin  6 cos  0 
Ky  = K sin  0 sill  0 
K = K cos  0 

Z 

k = Vk  ^ + K ^ + k 2 = Ik  I 

f X y z — 

The  magnitude  .sign  on  k is  meant  to  change  a vector  to  a scalar  and  not  to  lake  the 
magnitude  of  a complex  number.  Thus,  k can  still  be  complex.  We  will  now  intro- 
duce the  direction  cosines  of  the  k vector  (a, /3,  7),  which  can  he  determined  Irom  the 
equations  above.  | 
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If  I - the  III  llic  ilrirriiiiiiuiil  ol'  llic  ilirlciirir  IriiMir 


I*  I'J  {(y^^P^)  + e„(7^ 

+ |e^^a7+e^^J7|  t €^^Py^or*ay  |a(3e^^. +^7e^  J 

+ e^.^P^ya  + ay  |e^^(7^+a^ ) - ^ f a^)  07 

^ -e,,ta<31  * 0/3  (/3^ +a^ ) 

+ (y^+P^)  Py-e^^P^y^l 

A sdliilidii  to  ('i|iui|joii  (',i'2)  will  I'ivc 


K,  ,=  ](JI)) 


-1 


[11  ± /iT^ 


•I  k I)  lel 


(3:J) 


wlicn-  wf  iioU'  two  [xissililc  cluii»cs  of  sipi  and,  ln-iur,  associate  llu-  suliscripts  1 and  2. 
Actually  then*  arc  lour  solutions,  with  the  possihilitv  of  a pins  or  minus  sipi  in  front  of 
each  ol  the  K’s  >fivcii  hy  (23).  The  sipis  in  front  of  the  K’s  are  ineons(‘(pi<’ntial  lieeaiisi* 
they  only  determine  the  direction  of  propagation,  which  will  la-  calculated  in  a 
houndary  value  prohletn  hy  the  actual  houndary  conditions.  The  solutions  for  K | and 
K2  };iven  hy  (23)  are  implicit,  since  the  elements  of  the  dielectric  tensor  are  functions 
of  K.  In  order  to  oiitain  an  explicit  solution  for  Ki  and  K^,  we  viill  perturh  with 
respect  to  e^.  As  a first  a|iproximation,  let  e.^  equal  zero,  then  K]  and  Kj  are 


k/'>  = k 


K,/'>  = k 


(24) 

(25) 


riie  sujierscript  one  indicates  a first  approximation.  A second  approximation  can  he 
obtained  hy  usiii}'  (24)  and  (35)  in  the  dielectric  ten.sor  and  hy  lettiiijj  he  finite,  hut 
not  too  lar^e.  Then  the  second  approximation  ('an  he  ealculati'd  hv  using  (32).  In 
general,  the  mean  wave  in  the  random  medium  will  have  the  form 


<K  (rp>  = A,  e'j-'-+  Aj  e'^-’”- 


(26) 


K,  - K,  (a,  a^  + (3,  + 7,  a^) 

!^2  ~ 1^2  (®2  i*x  ^ ^^2  Sy  72  *5,  ) 
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W hell  eomputiii^  the  elements  of  the  dieleetrie  tensor,  the  direetion  cosines  of  each 
wave  will  he  computed  usin^  the  first  approximation  for  K,  and  K2  alonu  with  the 
appropriate  houndary  eonditions. 

Kaeh  wave  in  (36)  must  satisfy  eipiation  (31)  hy  itsidf.  This  fact  can  h<-  used  to 
determine  the  relationships  ainoni;  the  electric  field  components. 

''lx  = «lx  'iz  (■■<") 

'ly  = ^ly  'iz  m 

- ^12  i*23  - i*i3  a22 

djx  - 

aji  a22  - ai2  a2i 


- “21  ^13  - aji  a23 

- 

^22  i*ll  ■ ^12  **21 


Tlur  coefficients  a.^  are  the  elements  of  the  square  matrix  in  eipiation  (31).  with  k 
heini;  set  eipial  to  K| , a = a, , /J  = /J, , and  7 = 7, . A similar  relationship  can  he  written 
for  and  A^^, . 


'2x 

= ‘•2x 

'2z 

(3<>) 

'2y 

= ‘>2> 

'2z 

(40) 

- 1'12 

1*23  - 

1*13 

l’22 

I'll 

1)22  - 

1*12 

1)21 

- I'21 

^23  - 

»*.3 

1)22 

l>u 

1)22  - 

hi2 

bji 

The  coefficients  hj.  are  the  elements  of  the  sijuare  matrix  in  equation  (31)  with  K heiii}! 
set  equal  to  K2 . Tlie  general  equation  for  the  mean  wave  now  hecomes 

<E(r)>  = A,^  (a,^  9x  ^'^ly  ^ + 

''2z  ^^*2x  -X  '''  '*2y  -y  ^ ^ ' (+•  ) 


()iu‘  fould  now  us*!  tiu!  I'onu  *j1‘  the  nifaii  wav*-  •'ivfii  1»_\  (H)  ami  |»lai-«-  il  in  (20)  to 
(IftiTinin*'  llu‘  s(-att(T(‘(l  wave  in  a raniloin  incilinin  of  infinit*'  (‘xicni.  It  slionlil  U- 
clear,  huwi'vcr,  that  *‘qnation  (41)  was  derived  witliont  inakin;;  any  assuinplions  aliont 
till'  direction  of  propaj^ation  or  the  form  of  the  correlation  function.  \\*‘  al."!)  a.ssnmed 
that  tlu‘  thr*:e  eipiations  pven  liy  (31)  w*’re  coupled.  We  will  now  <io  hack  to  eipiation 

(31)  and  .S(‘e  what  happens  wlum  the  direction  of  propagation  of  the  mean  wave  and  j 

the  correlation  function  of  p (r)  take  on  spi-eifi*-  forms.  If  the  a.ssumption  of  de 
pendenee  is  correct,  one  can  .s«‘e  from  (41)  how  depolarization  ini'ihl  ari.se.  I'or 
instance,  suppose  a wave  that  is  polarizi'd  only  in  y is  ineiihmt  from  fn'e  spae<*  onto  a 
half  space  of  random  media.  Then,  from  (41)  the  mean  wave  would  venerate  x and  / 
components  as  well  as  a y component,  and  a depolarized  term  in  the  seattereil  field 
would  result. 

Now,  let  us  consider  what  happt^ns  to  equation  (31)  when  we  a.ssume  that  the 
correlation  function.  15  (r  - r'),  is  iscjtropie.  'I’lus  means  that  15  (r  - r’)  is  a funetion'oidy 
of  I r - r'  I and  is  therefore  ind*;penden'  of  direction.  W e will  also  a.s.sutne  that  the  direc- 
tion cosine  of  K with  respect  to  the  y axis  (j3)  is  ecjual  to  zero.  This  assumption  is 
certainly  allowable  since  we  can  orient  our  coordinate  system  in  any  mannt'r  we  wish; 
therefore,  we  will  choose  tlu  crientation  that  .s(;ts  p e(|ual  to  zero.  15y  iisiii};  the.se 
assumptions,  the  element  ot  the  dielectric  tensor  is 

= K.’  </ f,  .J,'  (42) 

T ' 

w here  K = K (a  a^  + 7 a^  ) 

Lettiiiv  X - x'  = u,  y - y'  = v,  and  lettiuj;;  z - z'  = w,  and  rememherinv  that  v'  is  over  all 
space,  equation  (42)  hecomes 

^xy  ^*^0^  f J 1^1  '’3  (d)e^-’“  (^;5) 

w h(!re  K = Ir  - l’  1 

The  uppi  r part  of  fifiure  3 shows  the  relationship  hetween  the  K and  R vectors  in  the 
u.  v.  w coordinate  .system.  Since  /3  = 0,  the  K vector  is  in  th«'  11.  w plane.  We  want 
to  rotate  the  u,  v,  w coordinate  system  to  u'  , v'  , w’  in  order  that  the  w'  axis  coin 
eides  with  the  direction  of  jx.  'I’he  resultant  rotated  coordinate  svstem  is  shown  in  the 
lower  portion  of  fivure  3. 
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Figure  3.  Geometry  of  K and  R Vectors 

The  niathcrnatical  relationship  between  the  unprimed  and  primed  eoordinates  are  pven 
Ixdow; 


n = u'  cos  0,^  + w'  sin  0^  = iT  7 + w ' a 
w = - u’  sin  6 + w’  e()s  0 = - 11'  a + w ’ 7 

Z Z ‘ 


K({uuti<>M  (4il)  ciiii  How  l)<‘  wrillfi)  as  lollows: 


.ao  .00  ~9o 

liv  cliuiijiinji  to  s[»lH'rical  tooriliiiatos,  llu*  above  et^ualion  becomes 

2 rr  ff  » 

e - ^ j (i<p j <10 j i\\{  H^siiiOIHH)|'^i"Ositi0(7siiiO«os0  + a(  osO)|l'3(K)<' 


jKKru.>/? 


u 00 


Carrying  out  llu'  integration  for  tlie  al)o\e  e(|iiation  in  0 shows  that  is  eijual  to 
zero.  A similar  analysis  will  show  that  is  also  zero.  I hese  results  will  make  •‘(illa- 
tion (31)  appear  as  given  helow  : 


ait  I ) a A 

“ 13  * 

t)  a22  h A - 0 


0 a„  A 


where 


an  - 7^  - k„^  e 


O XX 
. J.  I-  2 


ai3  = - (K^7«  + ko^  ^xy.) 


a22  - - k^^ 


Expression  (44)  shows  that  the  (equations  become  uncoupled  by  allowing  the  correla- 
tion function  to  be  isotropic.  This  means  that  the  equation  derived  for  the  mean  wave 
earlier  (41)  is  not  valid  for  an  isotropic  correlation  function  and  that  a new  mean  wave 
must  be  sought.  For  a wave  polarized  in  y such  that  A^  0,  we  must  have  ^22  - 
This  condition  yields  the  following  solution  for  K,  which  we  designate 


^yy 


K„  = 

E =1  H(K)|f2(K)+  f3(R)|cJ^'<i-- 


dr' 


¥ 


Tlif  dcfiiiilion  of  I.  ^ivcti  aliovc  coiin-s  from  c . I’lic  value  of  K lo  he  used  in  I,  when 

* y y — 

eoinputiiii;  ( fo)  would  lx-  I’lie  vailuine  inlefiral  is  evaluated  in  appendix  A alon*j 

with  the  other  elements  ol  the  dieleetrie  tensor  and  the  direetion  cosines  a and  7. 
It  can  also  he  s«'en  Irom  ('(piatioii  (45)  that  a first  approximation  for  K (which 
occurs  when  e is  e(pial  to  zero)  will  simply  he  k.  We  will  rewrite  K as  follows: 

y y 

^yy  ^ ('1  + jyi)’'' 

X,  = Ke  I I + e/  l./k^} 

y,  - Im  |l  + k^''  e/  L/k^} 
l$y  ehanijin};  lo  polar  coordinates  we  have 


•^yy  - k <'xp(j0i/2)  (46) 


01  = tan''  (y,/x, ) -7r<  0,  < tt 


K({uation  (46)  represents  an  expression  for  the  effective  propajjation  constant  when  the 
mean  wave  is  polarized  in  the  y direction. 

I 

II  we  allow  the  mean  wave  to  he  polarized  in  the  xz  plane  such  that  and  A^ 
are  nonzero,  then  the  following  determinant  must  be  set  equal  to  zero  in  order  lo 
obtain  a nontrivial  solution  for  K: 


(let 


an 

ai3 


ai3 
a 33 


_ 2 

~ a|,  a33  - ai3 


0 


The  above  equation  will  yield  the  following  scjlution  for  K.  which  we  designate  : 


K...  = k 


I e e - e ^ 

f XX  zz  xz 

i ^xx  + T'"  + 27a;e^,_ 


'/i 


(47) 


As  a first  approximation  for  we  set  equal  to  zero.  When  this  is  done,  the 
following  re.sull  is  obtained: 

K = k 


1 

i 

i 
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l{y  clianf^iiif!;  to  polar  coordinates,  we  have: 


02  = tan'*  (ya/xa)  -Tr<(j>2<Tr 


We  have  obtained  two  solutions  for  the  effective  propagation  constant.  One  solu- 
tion is  used  when  the  mean  wave  is  polarized  in  the  y direction.  The  other  solution  is 
u.s«‘d  when  the  mean  wave  is  polarized  in  the  xz  plane.  The  next  section  of  this  report 
will  consider  the  problem  of  a horizontally  polarized  wave  incident  obliquely  at  an 
angle  6-  from  a free  space  medium  (z  > 0)  unto  a lossy  random  medium  (z<0).  The 
boundary  conditions  for  the  mean  wave  will  be  satisfied;  then  the  scattered  field  in  the 
free  space  medium  will  he  calculated.  Finally,  an  equation  for  the  radar  hackscatter 
coefficient  will  be  derived.  Following  this,  a .similar  derivation  will  he  performed  for 
the  case  where  the  incident  wave  is  vertically  polarized. 

Horizontal  Polarization  Analysis.  The  results  of  the  previous  section  on  Scatter- 
ing Geometry  and  Wave  Propagation  In  a Random  Medium  of  Infinite  Extent,  will  now 
be  used  to  derive  a radar  hackscatter  coefficient  for  the  case  of  an  incident  wave  which 
is  horizontally  polarized,  ('.onsider  figure  4,  which  shows  a plane  wave  incident  from 
free  space  (z  > 0)  at  angle  0.  onto  a half  space  (z<  0)  of  random  media.  The  polariza- 
tion of  the  incident  wave  is  in  the  +y  direction.  The  coordinate  system  is  right-handed 
so  that  the  +y  direction  is  into  the  plane  of  the  paper. 


24 


Figure  4.  Geometry  for  the  Coherent  Waves  (Horizontal  Polarization) 

The  solutions  for  the  mean  waves  in  tite  two  rejrions  will  be  considered  first.  The 
scattered  waves  will  he  considered  later.  The  total  mean  electric  field  (E, ) in  the 
upp<;r  medium  can  he  written  as 

E,  = a^  +R^eAl!r-lj  (4«)) 

n.  = a sin  6.  - a cos  0 

-1  -X  1 -Z  I 

n = a sin  0.  + a cos  0. 

— r —X  I — z \ 

r = X a + >•  ii  + / a 
— X -y  ^y. 

is  a reflection  coefficienl. 


1'lie  first  term  in  (49)  n'presents  tlie  iricidt^nt  wave,  and  llie  second  term  re[)resents  tli(^ 
reflected  wave.  The  mean  electric  field  (£2)  in  the  lower  medium  can  he  written  as 


Fo  = a 

—2  -y 


-j  K n,  • r 
e * yy  - 


(.-,()) 


113  - itj  ■"‘in  4^  - a.,  'P 


j 
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(51) 


'1'^  is  a transmission  ttH^fficitml. 
i//  is  llic  i omf)l<‘x  an*;!*-  of  refrac  tion. 

From  Siu;ll’s  law,  the  following  c'xprcssion  holds: 

K sin  li/  = sin  0- 

yy^O  I 

Hy  following  Stratton^' , an  cejuation  can  b<‘  dc-rived  for  cos  \p.  If  we  write  ^ as 
j , then  sin  4/  is 

sin  i/'  = (a  + jh)  sin  6^  (51 

where  a = ) 

and  ij  = k„ay  (^j,^  +ay2) 


If  we  solve  now  for  cos  \p,  we  have 


cos  \p  = V 1 - (a  + jh)^siir  0^  = e (5;{) 

The  magnitude  Pj  and  the  phase  <p^  can  he  found  hy  sejuaring  (53)  and  ecpiating  real 
and  imaginary  parts  on  either  side  of  tlie  equation.  The  subscript  t lias  been  used  to 
indicate  reference  to  the  mean  wave,  which  is  transmitted  into  the  random  media. 

cos  2 0,  = I + (h^  - a^)  sin^0j  (54) 

p^  sin  2 0|  = 2 ah  sin^  0^  (55) 

Solving  equations  (54)  and  (55)  for  p^  and  0,,  we  have 

Pt  = I ( 1 + (h^-a^ ) sin^  I ^ + da^h^sin^dj } (56) 

( 2ah  sin^  0j  ) 

0,  = 14  tan"'  ■( — r T ^ ) (57) 

* I 1 + (h^-a^)  sm^0i  ) 

An  exprc:.s.sion  can  now  he  written  for  K cos  0 as  follows: 

r yy 

cos  0 = (|3^,  - jofy)  p,  (cos  0,-  jsin0,) 

cos  0 = (]  - jp  (5») 

julius  Adams  .Slralton,  F.lpctrftm/tfinoUc  Tbefiry*  MrUraw  Hill.  lO  U . pa^zr  .Adi . 
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wluTf 


P = P,  (Py  Si"0,  +«y  COS0,) 

M -Pt  (Py  -“y  ^t) 


(59) 


'I'lif  iiu'aii  electric  field  in  the  raiidoiii  inediuiii  can  now  he  written  as 

K.  (r)-  (61) 

Tlie  tangential  components  of  the  electric  and  niafpietic  fields  must  he  continuous 
across  the  boundary  at  z = 0.  The  majinetie  fields  in  each  nunlium  can  he  computed  hy 
using  the  electric  fields  atul  the  Maxwell  equation  V x K = - jco/i  Ij.  Tlu!  houndary 
conditions  will  permit  a solution  for  Kj^  and  T^. 

k„  cos  0 . - K cos  \b 

K ^ (62) 

cos  0.  + cos  \P 

2 k cos  0- 

T,  = r ‘ (63) 

^ K cos  li/  + k cos  0. 

yy  o 1 

We  are  now  in  a position  to  compute  the  scattered  fields  and  then  to  calculate  the 
radar  hackscattcr  coefficient  using  the  scattered  field  in  the  upper  medium  (free  space). 
Tlie  etpiation  for  tl>e  scattered  electric  field  in  the  random  medium  is 

££*=  K Ii> 

I V(  V'  4)  -V"  H - k^KJ  - k„^  fc^/i  (r)  a/i;  e-jV-‘“''®iei‘i'eP*  (64) 

If  we  say  that  is  of  order  in  magnitude,  then  V ' eiui  he  shown  to  he  of  order 

and  can  he  neglected  in  (64).  since  every  other  term  in  the  equation  is  of  order  e^. 
The  three  components  of  (64)  are 


<! 

K,  + k^  K = 0 

8X  SX 

(65) 

<S» 

0 

il 

€ (r)  T^  ,.-jl‘oXsm9i,.jqz,.pz 

(66) 

V 1 

h:  + k^  E = 0 

sz  s/. 

(67) 

.'solutions  for  and  can  he  written  in  Ihe  form  of  a Eourier  transform  as  follows: 


and 


-'K-'K 


•■'•«  < ^-y  '"•)  = (-0^)2  / /'*.  O'x ' >'> ) •■  ' »• 


'"x*  .>>>  dk  dk 

X > 


«dl) 


(<)'») 


In  tlic  al»(»vt‘  f(|iia(ion.',  |{^  and  1$^  an*  Ciinclions  of  tin*  Fonricr  variables  k^^  and  k^ . 
A solution  for  k ' can  easily  lx*  obtained  by  pnttinx  ((>*))  back  into  (67). 


k,  ' = Vk^  -k 


(70) 


Since  k is  complex,  eipiation  (70)  is  not  in  a fjood  form  for  ealenlations.  bettinii 
k = (3^  - . we  have 


k '^  = /3  ^ - 2j/3  a -a^-k^-k^  - p < 

L J^o  o o X y 


xJCV. 


(71) 


riie  rnapiitiide  p^  and  phase  0^  of  k,'^  can  be  found  by  .^quarinix  (7|)  and  e(piatin}i 
real  and  imaginary  parts  on  either  side  of  the  ecpiation. 


(72) 

(73) 


p - \/()3^-a;^-k^-k^)^+  4o;  ^ B ^ 

o X y ' o 

It  k^'  is  now  written  as  k^.  + jk.  where  k^  and  k^  are  both  real,  we  have 

k ' ^ k + j k. 

z r J I 

k^  = pj"  cos  (0^/2) 
kj  = p^  sin  (0^/2) 


■A  solution  for  is  much  njore  difficult  than  for  and  owinji  to  the  term 


(73) 


appearinj;  on  the  right  hand  side  of  ecpiation  (66).  1'his  term  is  a random  function  of 
all  three  coordinates  owing  to  p (r).  A general  form  for  the  solution  of  can  he 
written  as 
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(7.")) 


IMaciii^  «‘(|iiuti(>ii  (77))  into  (66)  uill  ^vt* 

-k^|{  + + k^  |{  i <‘\|)  (jk  X + jk  V)  <lk  (Ik 

(l>7r)^  / y I * > > > (1/.^  > ( ' ^ ^ ^ ^ 

-k  ^ e (x(r)  T,  (J**'  c'" 

O S i 

\\»!  will  now  (Icfiiie  the  Icrtii  S(kj^,k^, .z)  as  I’ollows; 


(76) 


’ - j «*xp  (-jl',' - jl^jV  ) (ly  (lx  (77) 


. . -jk  x?<ind- 
ju(r)t;  « ■ 


" (kyf  ‘••'y 


(7H) 


Tlif  assiiinptiun  lias  Ih'cii  made  that  (r)  is  Fourier  traiisforenahle.  l'siii‘j(7H)  in  (76) 
will  pve 

I 7 7i  'i'i5  I 

- 7 I / - k Ml  - k Ml  + + k"  II  + k M S (k  , k , z)  T,  e^'^"  eP"  x 

(27r)M  y I * y V y (Iz^  y o s v x’  y'  ^ i 


exp  gk^x  +jkyy)  dk^  dk^  = 0 


(79) 


One  way  to  he  certain  that  the  integral  on  the  left-hand  side  of  (79)  is  always  zero  is 
to  set  the  integrand  ecpaal  to  zero.  When  this  is  done,  the  following  second  order 
differential  equation  will  result: 


(Fll 

y 

ib} 


+ (k^  - k^^  - k^^)  11^  - - k„^  S(k^.k^,z)  ej-*^  e^^ 


(80) 


We  will  let  11^  he  the  sum  of  a complementary  solution  (Hy^)  and  of  a particular  .solu- 
tion (Oyp).  The  complementary  solution  can  he  found  hy  setting  the  right  side  of  (80) 
equal  to  zero. 


|{  =15  + |{ 

yc  yl  yz 


(81) 
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rii<‘  (|uuiititi(‘s  15^1  and  15^.,  arc  iii(lc|)cii(lctit  dI' /..  'I’lic  lorin  <»f  the  particular 
can  l)c  written  as 


liyp  = + vj 

I’lacin*;  (M2)  into  (Mt))  and  iisinu;  tlic  inctliod  of  variation  of  |>arainct»Ts  will  yiidd 


k T . , 

V-  S(k  .k  ./)  ..J'*'-  C*"  C-J^* 

•>;  L ' ' X • V ' 


•'ll 

dz  2j  k,^ 


i!'2  = 
(1/ 


Iv  ' e 1 


= S(k  ,k  . /)  ' 


(h::) 


(«d) 


(M4) 


lntc»;ratinjj  Moth  sides  of  (Md)  between  the  limits  of  a and  z will  >iive 

2 r X 

d?  + v,(a)  (85) 

a 

where  the  lower  limit  a is  a constant.  Integrating  equation  (84)  ladween  the  limits  of  h 
and  z will  jnve 


(86) 


where  h i.s  a constant. 


The  solution  for  can  now  he  written  as  follows: 


‘"•d’*'*'’”'  (2,) 


* a ""-J  z _ oo  J 

-oo  2 

e + l^v,  (h)  + I S(k^ , ky , ^)  e^''^  e e''‘  ^ /s(k , . k^, , ^)  e e x 

d^  je'^'^*''  I exp(jk^x+jk^,y)dk^  dky  (87) 


hh\' 


* + 8^2  e-j^'-^  + |v,(a) 
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wluTc  liavf  madf  the  iiitffiral  in  v,  (/.).  a sum  <»f  lw<»  inlftirals.  Siinilarlx.  \j  {■/.) 
lias  lii'fii  made  <‘<|ual  lu  the  sum  <d’  two  inlc<crals.  \\<-  will  now  define  the  eoeffieienls 
lii  and  ISj  as 


15,  = 15^1  +vi(a) 


/s(k  .k  .^)e^"^  e‘’‘  dj 

k > 

J £ a 

152  =15,2  +v,(l.)+  j S(k^.k^.?)e'''S-‘’^-'^'‘  dj 


* I. 


Ineorporatina  (he  aliove  definitions  into  eijuation  (157)  will  yield 

z 

+ -'  /|^y^S(k^.k^,,Of(^)  e'J‘‘^^|  exp  (jk^x  + jk^  y)  'Ik^dk^.  (HH) 


where  f(?)  = eJ**^  e^^ 


Tile  first  term  in  the  inte};rand  of  eipiation  (Md)  represents  a s<‘t  of  waves  propa<;alin>: 
in  the  minus  z direction.  The  second  term  represi'iits  waves  propa*£atin<:  in  the  jilus  ■/. 
direction.  From  the  radiation  condition,  F,^^  must  }£o  to  zero  as  z approaches  minus 
infinity.  'Phis  means  that  152  must  hi!  zero.  The  rest  of  the  terms  in  the  integrand  do 
f;o  to  zero  as  z approaches  minus  infinity,  although  the  last  term  does  so  only  hy  using 
l/llospital’s  rule.  The  final  expn^ssion  for  is  then 


K (x,y,z)  = 


(l^x’  •'y-  ^ I ''-'f’  j'^y>  ) '"'x'*''y 


(27r)^ 


2''r^  T 

//  j 15,  S - ^)‘(?) <u] 


(»')) 


The  scattered  electric  field  (K^')  in  the  upper  medium  (z  > 0)  must  satisfy  the  follow- 
ing homogeneous  wave  (‘(|ualion: 


;{| 


1 


v'  K 


t k 


I-;  ’ 
-»* 


I) 


vom|)t>iifnt>  of  K ’ cati  1m‘  written  as 


K ; (x.y./.)  = , I /\,'(k^,  k^ ) .-x,.  Ok,x  + jk^  y - jk,/.)  .Ik^.lk^  (<>()) 

(2^p  / > ('»t) 

oo  oo 

K (’‘'V./.)  = J j\  (k^.  k^ ) <-x,,  (jk,x  f jk^  y - jk,/.)  .Ik^.lk^  ('JL>) 


I luTc  are  six  unknowns  assovialetl  with  the  eotn[)onents  of  the  sealltrred  ekn  trie  fielil> 
iti  tlie  two  mediums.  I'hese  unknowns  are  A^.  A^,  A^.  If^.  If, . and  IJ^.  Tlierefore.  six 
independent  erpiations  must  la*  devtdoped  that  will  allow  for  the  solution  of  the  six 
unknowns.  Aelually,  we  are  only  interested  in  solutions  for  .\^.  A^.  atid  . whiili 
will  he  u.s«(|  to  eompule  the  far  zone  scattered  field.  At  z = I),  the  taiifiential  eom 
ponents  of  the  eleetrie  and  ma^nelie  .seattered  fields  must  he  eontinuous.  This  can  he 
slate<l  as 


n4.-iis  I = 

(‘H) 

where  jj  , and  lj_.  are  the  scatlereil  maf/tnetic  fields  iti  the  upper  medium  and  in  the 
random  medium,  respectively.  I'hc  two  houndarv  conditions  fiiven  above  can  he 
restated  in  terms  of  the  electric  field  eompotienis  as 


I'Jfx.y.O)  = 

''sx 

(x.y.O) 

(‘h>) 

lV(x.y.<>)  = 

■■'■sy  ' 

(x.y.tt) 

<<Ky) 

an  ' 

ai-:  ' 

aiv 

ar. 

s/ 

— 

3z 

dy 

- — at  / = (1 

3/, 

ai;^ ' 

d\\; 

di: 

ar. 

— l*- 

— zz 

s\ 

- ~ at  z = 0 

C'H) 

fl/ 

d\ 

0'/- 

ax 

\2 

In  udditioii  to  tho  f<nir  »‘i|ualioiis  ^iv«‘n  above,  there  are  two  diver{;»‘nee  eoiiditions  that 
ran  he  us«‘d.  The  diver«enee  of  the  scattered  electric  field  in  the  upper  niediuni  is  /a'.to. 

V k;  = 0 

k A +k  A -k  A =0  (99) 

XX  y y i z 

The  divergence  of  the  scattered  electric  field  in  the  random  medium  will  he  zero,  pro- 
vided we  use  only  the  complementary  solution  for  If  we  consider  only  the  com- 

plementary solution  for  then  the  scattered  fndil  iti  the  random  medium  obeys  a 
homogeneous  wave  equation  of  the  form 

K + k^  K = 0 

v —sc  — 8C 

The  subscript  c has  been  added  to  indicate  that  only  complementary  solutions  are  to 
he  considered.  The  above  wave  equation  implies  that 


k H +k  B,  + k,;  II  =0  (100) 

Evaluating  equations  (68)  and  (90)  at  z = 0 and  placing  the  results  in  the  boundary 
condition  given  by  (9.5)  will  yield 


1 

(27r)3 


/ / { 'x 


x+jk^yOdk^  dk^  = 0 


Taking  the  Fourier  transform  of  the  above  equation  will  give 

\ = '^x 

When  eajuation  (89)  for  (x,y,z)  and  etjuation  (91)  for  (x,y,z)  are  evaluated  at 
z = 0 and  placed  into  the  boundary  condition  given  by  (96),  the  following  ecpiation 
results: 


-B,  + 


0 ^ 

y^S(k^,k^.?)f(?)eJ‘‘*^  d?  >exp(jk^x-Hjk^y)dk^dk^  =0 
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rakiiif'  the  Fourier  Irunsl'oriii  of  Itie  al)ove  <-.\|)re.ssioii  v>ill  ijive 

k ^ e T, 

• A, -i«.  + (I. -I2)  = " 

“J  "'  z 


U 

wliere  i,  = J S k^ , J)  f (J)  e ^ ilf 
0 

aiul  \,  = I S(k^.k^.J)f(^)e''‘/^  .1^ 

. g0 

Siil»stitii(iii>i  flic  proper  e<piations  for  file  <*leetrie  fields  into  tlie  iioundary  eundition 
fliveii  l>v  (*)7j  ai\<l  evaluating  (lu-  final  result  at  ■/.  - 0 will  };ive  the  followinfi expression; 

— s I fl  jt^  A + ik  A -jk  B +jk'|{,  - I, 

• k ^ e , 'I'  ( 

'2  I ‘‘XpOl^x^  jl^yy)  ‘^*^x'*'^y  "" 

When  wf  take  the  Fourier  transform  of  the  equation  above,  the  followiiif;  result  is 
obtained: 

l^y  \ + l^z  S - k,.Hz  + ('.  + W 

Substitutinj;  the  proper  equations  for  the  electric  fields  into  the  houndary  condition 
wiven  by  (9B)  and  evaluating  the  result  at  z = 0 will  give  the  following  equation: 

■j*'--'  ^ 

.00  _oo 

Taking  the  P’ourier  transform  of  the  above  equation  will  give 

k A + k A + k ' U - k B =0  (104) 

F^quations  (99)  through  (104)  are  six  equations  with  six  unknowns.  The  quantities 
of  interest  are  A^^,  A^,  and  A^.  When  these  quantities  are  .solved  for,  the  following 
expressions  result: 
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^C'x’S)  = ^ ‘^x2(''x-K)'2 

wlurc  A^,(l'x'''y)  = - 


(lOr,) 


uiul 


wlier 


und 


MV-'x'-y. 


X . 

L'',' J 

L 

II 

X 

(k.,k^)l,  + A^2  (^^^)l2 

r 1 

\ ^^^  + 2k;(k^^  . k.k,')  1 

2j 

_k;(k/+k;+k>;)(k^  + k;)J 

A.(l^x’S)  = \l(l^x’S)'.  + \2(^-^)'2 

where  A^j(k^,  k^)  = [k^A^,(k^,  A^,  (k^,  ky)|/k^ 

and  A^2(k,,  k^)  = [k.A^^C^x’  S)  + S 'VC^x’  K >1/^ 


(106) 


(107) 


A solution  for  the  scattered  electric  field  (K/)  in  tlie  upper  medium  has  now  been 
obtained  in  terms  of  the  intefp'als  {|  and  I2,  and  also  iti  terms  of  the  inte^ais  in  k^  and 
k . In  order  to  calculate  the  scattered  far  field  for  the  case  of  hackscatter,  we  will  use 
the  Stratton-Chu  intcfp-al  as  modified  hy  Silver. The  Stratton-Chu  intcfiral  will  Ik* 
changed  to  a form  that  will  allow  the  re.sults  of  the  previous  derivation  to  be  utilized 
in  computing  the  far  field.  The  scattered  far  field  (E^f)  in  the  direction  defined  by  the 
unit  vector jij  can  be  stated  as 


E^f  = ^ 112  \j  |n  X E^'  - T)  112  X (n  X Mj)|  - ds  (108) 

where  r is  a position  vertor  pointing  from  the  origin  of  the  (x,y,z)  coordinate  system  to 
a surface  element  ds  on  the  surface;  z = 0;  K = -jk^,  exp  (-jk^  H^)/47r8|j;  tlic  distance 

Samuel  Silver,  Microwave  Antenna  Theory  anti  Design,  New  Y ork.  Me(»raw  Hill.  194 1 , p.  I f»l . 
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is  llif  lo  llic  I'icld  point  ulicrc  is  caliiilalid;  K^'  and  IX.'  arc  llir 

total  scattered  eleetrie  and  ina^netie  fields  evaluated  on  tlie  surface  z = t);  t?  is  tlie 
intrinsic  iin|>cdance  of  frci?  space  (uf)per  inediuni);  and  n is  a unit  vector  normal  to 
the  surlaee  z = I).  'I'lie  surlace  intejtral  in  (I Off)  is  over  the  illuminated  area.  For  the 
i-a.se  of  hackscatt«‘r,  1(2  is 

n,  = - a sin  0 + a cos  0 = - n. 

The  mu^nclic  field  components  of  IJJ  can  he  computed  usiiift  equations  (90  to  92) 
aUmjt  with  the  Maxwell  equation  V ^ **  done,  the  follow  - 

ing compommts  of  II  ' are  ohiainetl: 

oo  oo 

(2Tr)K~  ff ^‘'v 

•oo  _oe 

oo  oo 

• OO  .OO 

oo  oo 

The  unit  normal  (n)  to  the  surface  z = t)  is  simply  the  unit  vector  in  the  z direction, 
w Inch  is  a,^.  The  calculation  of  the  (piantities  tii  x (n  x EJ)  and  112  x (n2  x (n  x )) 
are  now  straifjht  forward. 

92  (l>  X E^')  = - (a^  E^^'  cos  0.  + a^  E^^'  cos  6.  + a^  F,^^'  sin  0.)  (112) 

L‘2  X (92  X (n  X 11/))  = a^ll,y'  cos^  0^- a^  11^^'  + ajl^^'  sin  0jCos  0.  (lid) 

Flaciti};  equations  (90  to  92)  and  equations  (109  to  111)  into  (112)  and  (113)  and  eval- 
uating the  result  at  z = 0 will  give  the  following  result  for  the  scattered  far  field  (E^j); 

fjff  !-*  I + +’^x\V'^ol  + 

S -<» 

+ ay  [AyCOS0j  +(k^.A^  *^z''y)/’^o  I ■^il/.l'^x*'*'  + 


•^x  I (j*^x’‘  '*y 


We  will  now  allow  llie  limits  on  tin’  integrals  in  x and  y to  Ije  from  minus  infinity  to 
plus  infinity.  1 his  will  he  a rea.sonahle  assumpticin  as  lonj;  as  tin*  dimensiotis  of  the 
aetual  illuminated  area  are  all  mueh  greater  than  the  correlation  distaiiee  of  the  ramhjm 
dielectric  fluctuations.  'I'lie  integrals  in  x and  y will  then  yield  Dirac  delta  functions. 


where  is  the  incident  wave  and  a * is  u.sed  to  denote  a complex  conjupite.  The 
hrackets  aroutid  Kjm  are  u.sed  to  indicate  the  com|)utation  of  a statistii-al 

avera»«;.  The  quantity  is  the  illuminated  area.  Dy  u.sing  the  computed  form  of  A 
as  jjivei)  hy  (106),  E„„  can  he  written  as 


F 


0 


-jk  K — " 

L-  - /*  , f , r ly  , ^ -iik  xsmo.  jq£  pt 

4^il{ I ' ‘ <*^^(r)‘-  ■•■  <• 


Ayi(koSiii0;,())e'^‘‘*^  + 


Ay2  (k„sinOj.O)/'-^} 


(119) 


where  the  definitioiis  of  1,,  Ij  and  S (k^  siiiO.,0.  have  Iteen  used  in  (1  19).  The 
eoniplex  eonjugate  of  (119)  ean  now  he  written  easily  as 


-2jk^e'‘"'‘“eos0j 

47rR„ 


OO  oo  Q 

y'dx-y'  dy'  y’df'M(r')e""‘«’‘'“‘"''>e-^'^''e'’'' 


I Ayi  (k^sin0j,O)  e'‘*  ^ + A*^  (k„sin0j,  0)  e ^ } 


(120) 


An  ex|)res.sion  for  the  average  of  K|(||  Km,  ean  now  he  written. 

oo  oo  oo  oo  O O 


. .*  -jk't  . ,*  jk  £ Jk’*£'  , jk't  -jk'*£'l 

A ,A  ,e  * e + A ,,A  , e ^ e '■  + A .,A  e '■  e '■  I 

y\yZ  vzvl  v2v2 


y2‘'yl 


y2'y2 


(121) 


Till  term  </i(r)/2(r')>  which  appears  in  the  above  integral,  was  defined  earlier  as  the 
correlation  function,  IK[-[’)-  0i  our  case,  we  shall  let  H(r-r')  he  isotropic,  that  is 
IKE'I')  - 0(  I r-r'  I ).  Also,  in  order  to  make  the  mathematics  managcahle,  we  will 
a.ssume  that  the  form  of  the  correlation  function  is  exponential. 

l{(R)  = e-'‘/“  |{=lr-r'l  (122) 


where  2 is  the  (correlation  disfanc('.  Thestc  assumptions  ahoul  llic  form  of  the  correla- 
tion function  cannot  he  justified  analytically,  since  no  one  seetns  to  really  know  the 
nature  (jf  the  correlation  function  for  actual  vegetation.  The  a.ssumptions  arc  made 
mereiy  to  simplify  calculations  so  that  an  analytical  restdl  can  he  ohtained  (ptickly; 
then  the  theory  ean  he  compared  to  actual  experimental  results.  The  assum|ition  of 
using  the  isotropic  correlation  function  given  above  ean  oidy  he  justified  on  the  basis 
of  comparing  the  final  re.sull  with  experimental  data.  Also,  any  use  of  an  anisotropic 
correlation  function  will  ncsult  in  very  great  mathematical  complexities.  The  following 
change  of  variables  will  be  made  to  facilitate  carrying  out  the  integration: 
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X - X =n 


y - y = Vy 


riiis  cliaiif'c  of  variabifs  v^ill  allow  < K||||  lo  lx*  written  as  fullows: 

ao  OD  OO  Q Q ^ 

k%os^0.  r r C c c ^ 


k cos^t?  C C C C C / 
J 'ly'J 


-2jk  T,  sui9  jq(t-£')  p(£  + t')  t . -jk,£  Jk'*£'  -jk^£  -jk*  {' 

'•  '■  ^ Sl  S2'‘  •'  * 

^ . ..  Jk,£  jk/£'  , jk'£  -jk,7£l 

" a^2  \i  *■  ••  + \2  \2‘‘  " ‘ (i2;i) 

riif  intcfirals  in  \'  and  y'  a|>|>('ar  to  he  meaniiifjless  heeans*’  the  limits  extend  from 
minus  infinity  to  pins  infinity.  However,  the.se  inte>;rals  actually  form  the  illuminated 
area,  since  we  want  to  have  the  averaije  haekseattered  power  equal  to  zero  outside  the 
illuminated  area.  We  can  now  write  equation  ( I2d)  in  the  followin>;  form; 


(k^co.-^s ) \ r r r r 

<i„i,Shi>=  " J ■'i.J'sJ  « / 'If  I'K-irW)'- 


, -2jk  n sin0- 
•’  o 'x  I 


-h  -li 


w<^t)  mH)  i . .*  -jk'?  )K*i  , -jk't 

1 \ \ ' + A 1 A (*  ''  V '' 

' V I > I y 1 y 2 

\^2  ‘ '■  ' ^v2  ^v2  ' *•  1 


The  integration  in  J and  will  he  performe<l  first  hy  usine  th*'  followin';  ehaiifie  of 
variables; 

%.= 


I’hen.  the  differentials  d^  d^'  will  heeonie 


•U  'It?.,  'It?.' 
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the  JacDbiaii  **  is  •'ivoii  as  l»)ll<)w>: 


Kif^ire  5 shows  tlif  area  of  iiitcfiralioti  for  rarli  of  llir  two  si  ts  of  varialilcs, 


IL 

^iL 

(-h,  -h/2) 


(0,  -h) 


Figure  5.  Areas  of  Integration 

If  we  now  re|)resent  the  entire  integrand  of  (124)  l)y  f(J,  £').  then  the  integrals  in  { 
and  will  become 


0 0 

/•  /' 


0 


j j «ir  f(^  5')- 


dr?'  / (Itj  ffTj'  + '/2  T?  , TJ  - '/2  T?  ) 


-h  -h  \,(\) 

where  Tj^,  (t}'J  and  (v'j  are  defined  below 

^21  when 

r?22  (^1)  = 2 (T?>  h)  when 


- h < T?'  < - h/2 

- h/2  < T/'  < 0 
-h  < T?'  < - h/2 


- •'/-  ^ % < <» 


riu-  fiini  tioii  f(Tj,,  + '/2  ’?/'  \ - '’**  ‘■“l’"’'y  "*"■'■'■  '^2 

is  tin-al*T  that,  the  eorrelation  distance  C.  'I’hen,  the  intenration  is  only  over  tl.e  sl.a.led 
area  in  lij!:ure  5,  which  would  he  a few  correlation  leiifiths  wide.  I'lius,  tl.e  intcfiration 


in  T?  ainl  Tj  , now  heeotnes 
' '£  *• 


O <)  O « 

-I.  -h 

Tlie  above  a|)|>roxiination  of  the  inte<;ral  .should  he  valid  for  the  ease  where  h > C 
(correlation  distance).  The  equation  for  <F.„„  F,*„  > as  {^iven  hy  (124)  now 
heeonies 


K-„>^  A /•,,  rd,, 

^'^'11 II '^'11  It  ^ 4^2  ||2  J 'x  I 'y  I tj  X y 

Jt'-z  iv  V*  A* 

‘‘  ( 'yl  'yl  ^ yl  y2 

+ \ v*  +A  \*  (12.'-V) 

'v2  'vl*  >2  y2 


y2  'yl' 


The  inlcf'ration  in  tj*  is  now  quite  sini|de,  and  when  it  is  carried  out,  the  following 
e(|uation  for<  F||,|  F*||,  > will  result; 


(1^0  ‘ 

<F  F*  > = ^ 


, -jV‘z-''z 


'y2'yl"3 


where  e, . cj,  Cy  and  are  defined  helow: 


-’I'-jJv.-jlv, 

I 

-1‘^jK^jC 

I 


/fi<-  vol.mu-  integral  in  (|26)  is  In  sl  .-valnaU-.l  if  a rl, «..«»•  is  n>a.h-  lo  spheri 

ra  (•()(, rdinates.  ^ ifTiire  6 sinnvs  lln- -Iri.- n lali.)iis|,ip.s  l/rtunrii  , tj  , 1,  and  (lie 

splu-rn  al  coordinales  |{.  x,  0.  Alsn,  lh<-  .omdalion  function  is  assumed  to  bdsotropic 


and  exponential. 


|{(r7, 


- f T)  + 77  + 77  “ 

• ^v-’?,)  = e * y ' 


/e  -K/t 

= .•  =IJ(K) 


figure  6.  Spherical  Coordinate  Geometry 
Kquation  ( 1 26)  can  now  be  written  as  follows; 


“ n 2it 

f f •lx/'  d0l{^sinxe'''^'e'^'’‘o*'™'^'^®i  . 


jqRco?x 


V yl‘3'  +\2\2''4' 


jKrosx(k,, 


(127) 
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The  iiilejirals  in  (127)  eaii  now  he  worked  out  analytically  and  are  in  a|)|iendi\  15. 
After  earryin;;  out  tlie  integration,  the  following  equation  for  Kf, > is 

obtained: 


^*■‘1111  *■'1111^  “ 


2k^e^\  cos^e. 


>l"y2’2  ^ 

I 2 


A.,.A:.,e. 


1 1+1,2  £2)2  1 1+1,2  £2  p 


^yl''y2‘^3  ^ ^y2'^y2*^4 


ll+h2£2]2  |l+h2£2]; 


(12H) 


where  the  quantities  h^,  b, , bj  and  bj  are  defined  below: 


1,  = 
o 

y4k^sin2  0.  + lq-(|, 

b,  = 

V4k^sin2  6j  + [q  -(k 

' -k;rv2i2 

V4k%in2  0.  + [q  + (l 

1*3  = 

y4k^sin20,  + Iq  + (f 

4-^i<’:)/2]^ 

I he  radar  backscatter  co»‘fficient  can  now  be  easily  ol>tained  from  the  definition  ))ro- 
vided  by  equation  (1  18). 


<111 


\ \ 

yl  yl*^l  *yr‘y2'2 


A.  .A 


' ( |l+b2£2|2 


^ '^y2  Sl^'3  ^ \2^y2  'A 


|l+b2£2]2  [l+b2£2l2 


(129) 


hquation  (129)  is  the  result  of  all  work  in  this  .section.  It  musi  be  remembered  that 
the  .solution  a.ssurned  an  isotro|,ic  ex|ionential  correlati(,n  function.  Also,  the  effective 
|iro[,af;ation  constant  for  the  mean  wave  was  derived  usin*;  the  solution  of  an  infinite 
space  problem.  In  order  to  be  strictly  correct,  the  eftective  propa^atitvn  constant 
should  be  obtained  from  a solution  of  the  half-space  problem. 


1 
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'I’ln‘  iifxl  sfctioM  of  this  rt‘|iorl  will  1m;  devoted  to  olituitiiii^  a solution  for  tin- 
haeksr-atter  eoeffieieiit  for  the  ease  where  the  iiieident  wave  is  vt'rtieally  polarizetl. 

Vertical  Polarization  Analysis.  I’he  results  of  the  section  <»n  S(  atteriii*;  (ieotm  try 
and  Wave  PropaMatiou  In  a l{aMdotn  Medium  of  Infinite  Kxtent  will  now  he  used  to 
ilerive  a radar  haekscalter  coefficient  for  the  case  of  an  incident  wave  which  is 
vertically  polarized.  We  will  now  consider  a wave  with  verti<‘al  jiolarization  that  is  inci- 
dent from  free  space  (z  > D)  onto  a his.sy  random  media  (z<  0).  The  {{eornetry  of  the 
situation  is  j£iven  in  fitinn-  7. 


Figure  7.  Geometry  for  the  Goherent  Waves  (Vertical  Polarization) 

The  |M)larizalion  vector  for  the  incident  wave  now  lies  in  the  xz  plane.  The  solutions 
for  the  mean  waves  in  holh  media  will  he  ohiained  first;  then  the  solution  for  the 
scattered  waves  will  l)e  obtained  the  same  way  as  in  the  horizrmlal  polarizatioii  case. 
Till-  total  mean  magnetic'  fiidd.  Iii(r).  in  the  upper  medium  can  he  written  as 


t-jk  n -r  -jk  ti  T 1 

0-.  ^ ..-r  - J 


n.  = a sin^?  - a vtisO 

- I - X I “ Z I 

n = a siii0  + a cos^ 

- r - X 1 • z I 

r = \ .1  V a + / a 

-X  - •>  -z 


1\  II  is  u rfllffliim  i-tH'l'l'icient. 

'I'lif  first  Icrni  in  (lift))  represents  the  incident  wave,  utid  llie  seeond  term  repn’senls 
the  reflected  wave.  'I'he  mean  ma^netie  field,  M2(r)'  ll*''  hiwer  medium  can  he 

written  as 

n3  ^a^^sin^  - a^^cutnj/ 

T||  is  a transmission  coefficient. 

\p  is  the  com|)lex  angle  of  refraction. 


\V  riting  Snell’s  law,  we  have 


k.  sini/'  = k sin0. 

XZ  ^ O 1 


Again,  hy  following  Stratton^^,  an  equation  ean  he,  derived  for  cosi/z.  If  we  write  as 

siin// = (a  + jh)  sinOj  (Idd) 

•'  = “xz  •'o^^xz  + “xz). 

Solving  now  for  eosi//.  we  have 


cosi// = y 1 - (a  + jh)^ sin^  6|  = Pj  e ^ ' (lifd) 

'I'lie  magnitude  p,  and  the  phase  can  he  found  hy  sejuaring  ( t d4)  and  equating  real 
and  imaginary  parts  on  either  side  of  the  equation. 

p^  cos20,  = 1 + (ti^  -a^)sin^0j  (lifo) 

p^  sin20,  - 2af)sin^  6-  ( 1 46) 

Solving  ('(piations  ( I 4.4)  and  ( 1 .46)  for  p^  and  0,.  we  have 

Pi  = -j  I I -(a^  - h^  )sin^  0|  I ^ ■l.a^t)^sin'*0|| (147) 


Jtiliij-'i  \(lain>  StraUon.  f-hu'tromafntf'tic  fhenry,  MrCiran  Hill.  New  York.  I*)4I.  pa^«’  501. 
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(^1  = '/a  lull*' 


2 a f)  0 
1 

1 - (a^  - 


All  r(|uatiiiii  can  now  lie  wrillcn  for  cos  \p  as  follows; 

K^^cos  tp  = /5,  (cos0,  - jsin0,) 

K^/os  i// = q - jfi  (|;{(>) 

where  p = p,  (a^^eos0,  + ( I Ml) 

ii'xl  h-  Pt  ) ( ' ' > ) 

The  mean  mafinetie  field  in  the  random  media  can  now  he  written  as 

,,  , , M/.  (h, 

ii2(l)  = ay  1||  ' e e (I  f_>) 

0 J_) 

h’rom  the  Maxwell  equation  Vx  ii  = J + - — . an  expression  can  he  ohtained  for  the 

o i 

mean  electric  field.  Kifr),  in  the  random  media. 

L’  / \ - * II  r I h sin0j  j -j>‘,.xsine.  pz  jqz 

"''  "'XT'  rl’Tri  T ‘ ’ 

The  transmission  and  reflection  coefficients  can  he  eoiiqiuted  from  the  tangential 
components  of  the  electric  and  majjiietic  fields,  which  must  Im-  continuous  across  the 
lioundary  at  z = 0. 

K„ , COS0.  - k eosii/ 

K,,  = -7^ ! ^ (IT4) 

K^Z  *‘Os0j  + cusip 
2K  cost?. 

O ^ XZ  I 

The  mean  electric  field  as  pven  liy  (143)  is  now  completely  determined  and  can  he 
used  to  calculate  the  .scattered  field  in  the  lower  mfdiurn. 


I V(V  t-.)- v'  ii'sl  = - 


K,  ' ll 


/ k >iu0  \ 

“> 

' \ y. 


-JO, 


-jk  p/  jqz 

,•  ■’  ‘ <•  .• 


(1 16) 


\\t‘  shall  iiffilfct  the  iirailk'iil  Ifrin  in  ( 116)  for  (he  same  reason  i(  was  in‘*;leete(l  tor  (he 
horizontal  [volarization  solnlion.  The  three  eom|>onents  of  ( 1 16)  then  heeome 


V 


,,  ,2.,  ''o  ' ll  -jo,  -jk^xsiiiOi  iV 

K^+k^  = r P,  '•  ‘‘ 


^ sy  sy 

kl  e^^X(r)u>^l^  1,1  -jk  xsinOj  ^i.  jq/. 

h -I-  k^  K = sm0.  e " ' e e 

V S'/.  SZ  Li  I 


(IlH) 

(146) 


\ solution  for  K ran  l)0  written  in  the  form  of  a Fourier  transform. 

sy 

ti'.-S)***'  '"-.'IS 


(150) 


In  the  above  equation  If^  is  a function  of  the  Fourier  variables  k^^  and  k^ . I'he 
(juantity  k^  has  the  same  definition  that  it  had  earlier,  as  fiiven  by  equation  (7()).  The 
solutions  for  F.  and  F are  not  as  straisihtforward  as  F.„  heeause  of  the  terms  on  the 
ri^ht-liaiul  side  of  (147)  and  (149).  These  two  terms  contain  /i(r),  which  is  a random 
futietion  of  all  three  coordinates.  By  considering  first,  a general  stjiution  can  he 
written  as 


F_,^(x.y.z)  - 


(27r)^ 


(151) 


Placing  equation  (131)  int(j  (117),  we  have 

d^B' 


OO  oo 


— —f  /*  I -k^  b'  - k^  b'  + — rr^  ^ "1  ( ‘■’‘I’ ^ i''v  y) 

(27r)^  J J \ ' ''z^  * ) ^ ^ ' ■ 


''o  'll  A -j^,  ,/l'-  „P* 

— P.  »‘ 


(152) 


XZ 
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II  we  n-call  llu-  (Ifliiiilioii  ol  S (k^ , . z)  as  hy  ( (jualion  (77).  llic  aliovf  < x|)n's- 

sion  i-aii  Itc  writfcii  as 


exii(jk^x+jk^y),lk^.lky  =0 


k^  lC-M„S(k,.k^.z). X 


wluTf  M = 


''  ll  *■ 


I lu  iiitefiral  on  the  lolt-liatul  side  of  (15.5)  van  bv  satisiivil  l)y  s/*ltin}'  Ihv  intv^and 
equal  to  zero.  When  this  is  done,  the  following  second  order  differential  equation 


results: 


^ Mk^-k^k;)H:.M„S(k„k^,^)e^^^^  (.54) 

I he  above  differential  equation  has  exactly  the  .same  form  as  equation  (80),  which 
was  .solved  [ireviously  in  eonnection  with  the  horizontal  polarization  case.  Applying 
the  same  variation  of  parameters  method  to  (154)  will  yield  the  following  solution  for 
'•'sx  <'-y’Z): 


-00  .oo 

exp  (jk^x  + jk^.y)  'Ik^dk^  (j55) 

where  Rj  is  a function  ol  and  k^.  Since  equation  (149)  for  has  exactly  the 
same  form  as  ( 147),  a solution  lor  E , can  be  written  as 


..  Pt  jqf 

t)  *’  e e * df 


where  N„  = 


•^o  ‘'ll 


XZ 


The  components  of  the  scattered  electric  field  in  the  upper  medium  (z  > 0)  can  he 
written  in  the  same  manner  as  was  done  for  horizontal  polarization  (e.*;.  equations 
90  to  92). 


(27r)^ 


If 


f^y(l^x'  ‘■’‘I’  Ol^x  ' + j^y  ■ '“‘x'^^y 


(157) 


(158) 


where 


2 

X -y 


Again  we  have  six  unknowns,  C,^,  Bj  and  B^-  I he  three  unknowns  of 

interest  are  and  C^,  which  determine  the  scattered  electric  fields  in  the  upper 

medium.  The  boundary  conditions  at  z = 0 that  will  permit  the  calculation  of  the 
wanted  quantities  are  as  follows; 


3E 

dE 

9E' 

9E' 

SZ 

^ 

— SZ 

.sy 

dy 

3z 

9y 

3z 

9E 

aE 

9E' 

dE' 

sx 

SZ 

_ SX 

3z 

d\ 

3z 

3x 

Ksy(x,y,t))  = K,y(x,y,0) 


at  z = 0 


at  z = 0 


(160) 

(161) 

(162) 

(16.}) 


In  addition  to  the  four  boundary  conditions  given  above,  we  have  Ihe  two  divergence 
conditions  that  also  can  be  ii.sed.  When  the  appropriate  equations  for  the  scattered 
electric  fields  are  substituted  into  the  aliove  boundary  conditions  along  with  the  diver- 
gence relations,  the  result  will  yield  six  algebraic  equations  and  six  unknowns.  The  six 
algebraic  equations  are  obtained  in  exactly  the  same  manner  as  they  were  in  the  hori- 
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zoiitui  polurizulioii  solution  anil  are  ^ivcn  liclovv. 


k C + k (;  =k  C 
XX  y y X z 


k It'  + k It'  + k'  It'  = 0 

XI  y y z 2 


(i.-i2)  = (:. 


It  = c 

y y 

k C +k  C -k  It',  +k'lt'  = (1,  -1,) 

y z z y V 2 z y oj,,-  ' > 

, , , (Ii  - *2)  M (li  + li) 

-j  •'z  “J 


when-  I,  =y' S (k^.  k^  . l (f)c 

O 

U=j  <1^ 


A . .A. 


W luTi  liquations  ( 164)  to  ( 166)  arc  solved  for  and  we  have 

I. 1= 


-k  k k M F,f 

V X Z <)  * ) 


C , = — I - k k k'  N F,  - 2k  k'^  M F, 

x2  *>11  K ' X z z o * z z o ‘ 

“j  /*  o 


-k^k  M F.  + k k k'N  Fj  + k k k M F,  I 

X z o 1 y z z o ^ y X z o ^ ( 


(171) 


c = c .1,  +(:  I, 

y yl  ‘ y2  ^ 


1 


2ik'F 

J z o 


j k k k'  N Fy  + k k k M Fa  - k k k'  N F^ 

(yZZO'*  yxzo-*  XZZO^ 


-kJk.M.F,) 

= iTTT;^  I .k,k.k.M„K,  +k,k.k;N„P. 

J *^Z  o 


+ 2k„k:"\LF,  +k"^k^M^Fa  \ 


Z Z o 


I,  I, 

i:.,  =(k,(;,,+k,c,,)/k. 


where  F =(k’k'  t k‘k'  +k’k  r k'’k  ) (k’k’  + k’k’  + k'’k  +k’k  ) 

o 'yz  zz  vZ  zz^'zz  xz  zz  xz^ 


-k"k^(k  +k')^ 

X y ' z z ' 

F,  =k^  k'  + k^k'  + k^k  + k'^k 

* y z z z y z z z 


(*2  - ^ 


Fa  = k\'  + k^k'  + k'^k  + k^k 

•’  Z Z X Z Z Z X z 


(172) 


An  I’cjuatioii  for  llic  far  zone  scattered  field.  F^p  was  developed  earlier  from  the 
Stralton-(dni  inlei;ral  and  has  the  following  form: 


tv  = 


-jk  It 

2jk^eos0je  " « 
AttH 


[a/:,(koSin0i,  0)  + ayCy(k„sinf?p  0) 


a/;,,(k„sin0p())] 


(173) 


'I'he  interest  is  now  in  the  x and  z eornponenls  of  ( I 73),  which  will  Torn)  the  vertically 
polarized  received  wave,  The  doul)le  subscript  is  used  here  to  indicate  the  polari- 
zation of  the  transmitted  and  receivcil  waves,  respectively. 
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K = K • (a  i'osd  + a 

vv  — S '—X  I -'I  l' 


2jk  cosO  f “ •’ 

K = fcosOYA  , I,  + A , L) 


+ sinO.  (A^,  Ij)] 

riif  above  fcjuaJidii  t an  l/c  mani()ulaled  easily  into  a form  that  will  allow  a (|ui(’k  solu- 
tion ft)r  the  baekseatter  eoeffieicnt,  a“^. 

2jk  eosti.e''‘‘"'^“ 

— 4;^ir (A,.  (k„smd,  (.)!.+ 


A,,  (k„sin0.,())l,|  (174) 

where  A^,,  (k^sin0j,())  = Aj^,(;osOj  + A^^  sin0| 

Av2  (k„sin0.,O)  = A^.cosO.  + A^^sine^ 

It  can  be  seen  that  the  form  of  as  <iiven  by  ( I 74)  is  the  same  as  that  of  Fm,  >iiven 
by  equation  (117).  Theref<>re,  a final  result  for  the  baekseatter  coefficient  0°^,  can  be 
written  directly,  since  it  will  have  the  same  form  as  ^°lll  ■ 


47rK, 


A„ 


< ^-;v> 


where  Kj  * £*  - (k^/toe^)^ 
a°  =8cj^e^  cos^0. 


Sa„a:,c,  ^ a,,a:,c2 

)li+b^e^p  n+b^cM" 


^ A^jA^jCj  ^ Aj,2Ay2’4 


ll+b^g^p  [l+b^g^l' 


(175) 


The  quantities  c,,  C2,  Cj,  (4,  b^,  b,,  1)2,  ba  in  (175)  are  defined  in  the  .same  manner 
as  they  were  in  the  Horizontal  Polarization  Analysis  section,  the  difference  Lkmii}:  that 
()  and  q must  be  replac'ed  by  f)  and  q respectively. 
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Rquatiori  (175)  is  the  end  result  of  all  work  in  this  section.  The  next  section  of  the 
report  will  show  some  calculations  of  a”,,,,  and  and  will  indicate  the  effects  of 
chanpng  various  parameters. 


DISCU-SSION  OF  IlKSULTS 

This  section  will  show  some  numerical  calculations  for  the  theory  derived  in  the 
|trevious  sections  and  will  indicate  the  effect  of  various  parameter  changes  on  the  hack- 
scatter  coefficient.  A com|)uter  program  was  written  for  the  solution  of  ecpiations 
(129)  and  (175).  There  are  five  input  parameters  to  this  computer  program: 

1 . Frequency. 

2.  Average  relative  dielectric  of  the  random  medium. 

.3.  Average  conductivity  of  the  random  medium. 

4.  Correlation  length. 

5.  Standard  deviation  of  the  dielectric  fluctuations  of  the  random  medium. 
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A «‘i)|)y  of  tlif  compulfr  |)ro>;rani  is  wivt-n  in  a|i|)<'n(iix  (1.  'I'lic  oiilpul  of  the  coin- 
|iut(-r  pro^ruin  is  the  hueksealter  eoeffieii'iil  in  (le('il)els  as  a fiinetion  of  the  incidence 
an>;le  for  the  two  (lolarizalions.  'I’lic  hackscallcr  cocffici«-nl  in  decihels  is  related  to  the 
hai'kscatter  coefficient  as  follows: 

a“  (in  decihels)  = 10  lo>;,o  a" 

The  hackscatter  eoeffieient  on  the  ri^ht-hand  side  of  the  above  equation  is  corn|iuted 
from  eipiations  ( 1 20)  and  ( 1 75). 

I 

Fipire  8 shows  a study  of  k^li  variations  for  horizontal  |iularization.  In  this  ! 

fifiure,  the  value  of  k^£  varies  from  0.5  to  8.0  with  no  average  loss  (o^=0).  It  can  he  j 

seen  that  as  k^C  increa.ses,  the  level  of  the  curve  goes  down,  hut  the  shape  remains  the  | 

same.  There  is  a large  drop  from  k^  8 = 1 to  k^6  =2,  much  larger  than  the  ehange  from  1 

k 8=2  tok  8=8. 

O O 

Figure  0 is  a study  of  k^8  variations  for  the  ca.se  where  kjj8  is  less  than  1 and 
a^=0.  .\s  k^8  increase's  from  0.2  to  0.4,  the  level  of  the  curve  goes  down,  hut  the 
.shape  stays  essentially  the  same.  We  .see  also  that  the  separation  between  the  curves 
is  approximately  the  same. 

Figure  10  pre.sents  a study  of  k^8  variations  for  the  case  where  has  the  finite 
value  of  0.01.  It  can  he  seen  that  as  k^8  increases  from  0.2  to  0.6  the  level  of  the 
curves  increases  instead  of  decreases  as  they  did  in  figure  9.  The  spacing  between 
curves  changes  also;  it  decreases  as  k^8  increases.  When  k^8  becomes  equal  to  I,  the 
level  of  the  curve  is  now  lower  than  the  curve  ‘‘or  k^6  =0.6.  Thus,  if  can  be  seen  that 
by  comparing  figures  9 and  10,  the  introduction  of  average  loss  will  greatly  affect  the 
final  behavior  of  the  curves. 

Figure  1 1 provides  curves  of  hackscatter  coefficient  versus  k^8  for  three  different 
values  of  with  the  angle  of  incidence  equal  to  10  degrees.  It  can  be  seen  that  when 
0^=0,  the  curve  is  monotonie  and  decreasing.  When  a = 0.001  and  0.01,  the  curve  starts 
at  some  low  value,  increases  to  a maximum  value,  and  then  falls  off  gradually  as  kjj8 
increases  further.  This  is  a very  radical  change  from  the  case  where  a^- 0.  It  is  not 
known  whether  this  change  is  valid  or  if  it  is  a consequence  of  the  assumptions  used  to 
develop  the  mathematical  model. 

Figure  12  presents  curves  of  hackscatter  coefficient  versus  angle  of  incidence  for 
five  different  values  of  the  average  dielectric  (e^).  It  can  be  seen  that  as  increases, 
the  mean  level  of  the  curve  goes  down.  Also,  as  increases,  the  curves  drop  off  faster 
as  the  angle  of  incidence  increases.  The  spacing  between  curves  decreases  as  e 
increasT's,  with  the  largest  difference  occurring  between  the  = 1 curve  and  the 
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Figures.  Study  of  Variations 
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BACKSCAHER  COEFFICIENT  IN  db 


POLARIZATION:  HORIZONTAL 


f = 9GHz 


Figure  10.  Study  of  Variations 
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POLARIZATION:  HORIZONTAL 


0.1  0.3  0.5  0.7  0.9  1.1  1.3  1.5 

VALUE  OF 


Figure  1 1 . Study  of  Variations 
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= 2 t Wlwii  iiurfascs,  llu‘  itia^iiitmlr  of  llic  n-flrcli-d  ciilii-mil  iiicaii  vka\r 

increases,  and  less  fiUTfiy  is  availalilc  for  l)arks«'allfriii>;.  I'liis  would  a|i|i<'ar  l<»  arroiaid 
for  the  ra|iid  di-cra-asi-  in  llic  liat-kscallfr  cordficifnl  as  itl(•rt■as(■^. 

l''i<:urf  Id  |irfsfnls  curves  of  liackscallcr  (■(u'fficicnl  versus  anfik-  of  ineidenee 
for  four  different  values  of  6 . tlie  standard  deviation  of  tlie  dielectric  fluctuation'. 

S 

It  can  lie  .seen  that  as  e decreases,  the  level  vif  the  curves  iroes  down,  and  al.'o  the  eurve> 
fall  off  faster  with  ineidenee  ancle.  The  distance  hetwfen  sueee.ssiva-  curves  d«-erea.se.' 
as  inerea.s<‘s  froni  0.1  to  0.4. 


Kicure  14  |)re.sents  curves  of  haekseatter  eoeffieietit  versus  ini  ideiiee  ancle  for  six 
different  values  of  o , the  tnean  value  of  eonduetivity  in  the  random  metlium.  As  o 

«1  it 

inerea.ses.  the  mean  level  of  each  curve  deereasr‘s,  ami  the  curves  fall  off  faster  with 
ineidenee  ancle.  The  S|)aeing  helween  the  curves  varies  as  in<-rea.s<-s.  with  the  larg<‘.s| 
sjiacinc  oeeurrinc  hi'tween  the  = 0 curve  and  the  o.^~  0.01  curve. 

Kicure  15  presents  curves  of  haekseatter  eoeffieient  versus  angle  of  ineidenee  for 
four  different  values  of  k^  C.  The  polarization  con.sidered  in  thesr*  curves  was  vertical. 
It  can  he  s«;en  that  as  k^C  inerea.ses,  tlie  mean  level  of  the  ••urve  g(H‘.s  down.  Also,  for 
.small  angles  of  ineidenee,  the  curves  increase  instead  of  deerea.se.  as  they  did  in  the  I'a.si- 
of  horizontal  polarization.  In  comparing  figure  1,5  to  figure  8,  we  sr-e  that  for  each 
value  of  k^8  .the  vertical  polarization  curve  has  a higher  average  level.  The  curves  of 
figure  15  do  not  drop  off  nearly  as  fast  as  the  curves  in  figure  8. 

One  of  the  primary  difficulties  that  has  not  been  addres.sed  is  relating  the  statis- 
tical parameters  of  the  random  medium  in  the  theory  to  the  actual  |)aramelers  of  the 
vegetation.  This  can  only  he  done  in  an  approximate  manner  at  pre.sent.and  we  shall 
approach  this  prohlem  hy  using  the  ecjuations  that  were  developed  in  the  introduction. 


Kquation  (5)  gives  the  value  of  the  relative  mean  dielectric  constant  and  is 
restated  Itelow: 


where  = the  volume  of  one  leaf  = w x w x I 

Nj  = the  total  number  leaves 

= relative  complex  dieleidrie  constant  of  one  leaf 
= volume  of  the  forest  which  is  solely  air 
= relative  dielectric  constant  of  air  = I .(I 
V.p  = the  total  volume  of  the  forest  = I,  x h x I) 

O)  = radian  frequency 

= permittivity  of  free  space 
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ANGLE  OF  INCIDENCE  IN  DEGREES 


Figure  13.  Study  of  Variations 
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BACKSCAHER  COEFFICIENT  IN  db 


POLARIZATION:  HORIZONTAL 
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Figure  14.  Study  of  Variations 
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Figure  15.  Study  of  Variations 


l ilt'  lor  (lit'  varianci*,  is  a>  lolluws: 


2 

e = 


+ V.v(e.^  -e/ 

V... 


wluTc  ej  is  llif  real  part  of  e^. 

Olio  way  in  wtiioli  could  lio  dotorminod  would  lio  to  use  tlio  rnodol  aivoii  liy  I’oak 
and  Oliver.^'’ 


wlioro 


!’■ 

Ko 

lin 


=(172)  Ko(e^,)+j(l7;{)  Irn(e^) 

tlio  fraction  of  water  by  wcifjlit  in  tlio  vojjotation 
tlio  relative  eoni|)le,\  dieleelrie  constant  of  water 
means  take  the  real  part  of 
means  take  tlie  imaginary  part  of 


(aisgrill,  Peake,  and  'I’aylor^®  presented  a model  for  the  relative  complex  dielectric 
constant  of  water,  which  is  a function  of  wavelensith. 


= 5 + 


75 


where  X is  the  wavelength  in  centimeters. 

If  we  now  chose  a value  for  F.  a wavelength  X,  and  the  size  and  separation  of  the 
leaves,  we  can  compute  values  for  e^,  a^,  and  e^.  It  remains  now  to  eom|)Ute  the  corre- 
lation length  C in  terms  of  the  dimensions  of  a leaf.  'I’his  was  done  by  Fung,^*  and  our 
derivation  will  follow  his.  riie  correlation  length  C is  a measure  of  the  average  size 
of  the  .scatterers  (leaves)  in  the  vegetation.  Since  we  assumed  the  correlation  function 
to  be  i.sotropie,  we  need  to  relate  the  volume  of  an  equivalent  sphere  with  radius  K.  (o 
the  volume  of  a leaf.  This  is  done  as  follows: 


i TT  = V, 
d * 


W X W X 1 
1/3 
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riif  corrchilioii  ilislaiicf  i)f  a .s|)liiTt'  can  l)c  shown  to  he  e(|ual  to  ().().'{  litnes  its  radin.'-. 
Therefore,  an  estimate  of  C would  he 


e ^ t).0;(  K 

S 


We  can  now  eoin|)nte  llie  haekseatter  coefficient  as  a funetion  of  moisture  content  in 
tlie  leaves.  As  a jiarticular  example,  let  us  consider  the  followin}!  parameters: 


h 

I) 

w 

t 

d 

f 


the  length  and  w idth  of  the  ve{;etated  re>;ion  = 15  meters 
the  de|)lh  of  the  vegetated  region  = d meters 

3 centimeters 
0.05  centimeters 
1 0 centimeters 

4 centimeters 

1.0 

frequency  = 0 (iHz 


The  total  numher  of  leaves  in  the  volume  can  be  calculated  as  follows: 


L 

to 

1) 

w + d 

a + t 

Using  the  above  e(]uations,  it  is  now  a simple  matter  to  calculate  , e^,  e^,  6^  and 

C for  any  value  of  F that  we  choo.se.  Figure  16  presents  curves  of  haekseatter  coeffi- 
cient versus  moisture  content  (F)  for  two  different  angles  of  incidence.  It  can  be  .seen 
in  both  cases  that  as  F increases,  the  haekseatter  coefficient  also  increases.  However, 
the  curves  increase  the  fastest  for  values  of  moisture  content  between  10  and  30  per- 
cent. As  the  moisture  content  increases  beyond  30  percent,  the  curves  start  to  level  off. 

Figures  17  and  18  provide  a compari.son  of  the  developed  theory  with  some 
experimental  measurements.  The  ex[)erimental  data  was  taken  by  Ulaby,^’  using  a 
radar  s|)ectrometer  that  recorded  measurements  from  8 to  18  (HIz.  The  particular  data 
used  in  figures  17  and  18  were  taken  at  a frequency  of  9 (Hlz  over  soyheans.  '1  he 
statistical  parameters  which  are  needed  as  input  to  the  theory,  were  siTiq)ly  estimated 
because  we  did  not  know  the  moisture  content  in  the  leaves  or  the  spacing  of  tlie 
leaves.  It  is  not  rigorous  to  compare  theory  and  experiment  in  this  manner.  However, 
we  can  obtain  some  qualitative  indications  about  agreement.  The  curve  in  figure  17 
compares  theory  with  experiment  for  horizontal  polarization  and  with  the  soybeans 
at  a height  of  56  centimeters.  For  angles  of  incidence  greater  lhan  20°,  we  have 
reasonable  (|ualitative  agreement.  For  angles  of  incident  e -s  lhan  20°, 


27  K.  T.  Ulaby,  T.  K.  Bi'rsh.  I’.  F.  Itallivala,  and  J.  Cihlar,  Radar  Responses  to  Vegetation  II:  S IR  (lllz  Band,  Uiiivcr- 
sity  of  Kansa<s,  l^awrence,  Kansas,  KSL  Technical  Heport  177-51,  July  1974. 
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BACKSCAHER  COEFFICIENT  IN  db 


POLARIZATION;  HORIZONTAL 


f=9GHz 


10  20  30  40  50  60  70  80 

PERCENT  MOISTURE  IN  THE  DECIDUOUS  LEAVES  (BY  WEIGHT) 


Figure  16.  Study  of  Moisture  Variations 
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BACKSCAHER  COEFFICIENT  IN  db 


POLARIZATION:  HORIZONTAL 


f =9GHz 


oa  = 0.001 
4oi=  1.2 

£s  = 0.2 


0 10  20  30  40  50  60  70 

ANGLE  OF  INCIDENCE  IN  DEGREES 


Figure  17.  Comparison  of  Theory  with  Experiment  for  Horizontal  Polarization 
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BACKSCAHER  COEFFICIENT  IN  db 


r ^ 


-15 

-20 

-25 

-30 

( 

ANGLE  OF  INCIDENCE  IN  DEGREES  I 

Figure  18.  Comparison  of  Theory  with  Experiment  for  Vertical  Polarization  j 
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the  experimental  data  is  considerably  higher  than  the  theory.  One  explanation  lor 
thi.s  is  that  al  small  angles  t)l  incidence  there  may  he  eon.sideruhle  ground  M atlning 
that  the  theoretical  model  does  not  lake  into  account.  As  the  angle  of  im  ideme  in 
creases,  the  ground  elTeets  .■should  diminish  since  the  ray  path  is  almost  totally  m the 


\egelation. 


We  conclude  thi.s  .section  with  a hriel  diseu.ssion  on  the  limitations  of  tlie  theory 
that  has  lieen  developed  in  this  report.  The  idea  of  simulating  a region  of  vegetation 
with  a eonlinuous  random  medium  .seems  to  he  valid,  hut  it  is  not  known  how  well  the 
renormalization  leelnhipie  solves  the  prohlern.  One  prohlem  is  that  we  don’t  know 
how  mueh  of  the  multiple  seatiering  effect  we  are  really  con.sidering.  .Anotlier  prohlem 
is  that  the  mean  wave  was  derived  Irom  the  solution  of  an  infinite  space  prohlem  and 
then  applied  to  a half-spaee  situation.  Kupiee,  ct  have  .shown  for  the  scalar 

prolilem  that  llie  mean  wave  in  the  pre.sence  of  a hounded  medium  cannot  he  a.scer- 
tained  merely  hy  a.ssigning  to  tin-  hounded  medium  the  effective  mean  wave  propaga- 
tion constant  for  an  unhounded  medium.  In  addition,  the  effect  of  a .small  transition 
layer  near  the  houndary  must  he  considered.  A solution  for  the  tran.sition  layer  effect 
for  the  vector  prohlem  does  not  exist,  .so  its  overall  influenee  on  the  final  solution  for 
the  haek.seatter  coefficient  is  unknown. 

Another  prolilem  with  the  pre.sent  theory  is  that  an  isotropic  correlation  function 
was  u.sed  to  reduce  the  mathematical  complexity ; hence,  it  was  found  that  no  depolar- 
ization results  were  ohtained.  Although  the  correlation  functions  for  real  \i-getation 
are  not  known,  they  are  prohahly  not  isotropic.  If  anisotropic  correlation  functioiis 
are  to  he  considered  in  the  solution,  the  resulting  mathmnatics  may  hecome  prohihi- 
tively  dillieult.  Also  throughout  our  derivation,  we  have  consistently  a.ssumed  the 
conductivity  in  the  random  medium  was  a constant  equal  to  o^.  Obviously,  this  is  not 
true  since  the  conductivity  changes  significantly  from  air  to  vegetation.  W hat  is  needed 
here  is  a random  medium  in  which  hoth  the  dielectric  constant  and  the  conductivity 
are  allowed  to  lluctuale  in  a random  maimer.  Once  again  the  mathematics  of  this 
more  realistic  ca.se  will  he  very  difficult  and  may  even  make  an  approximate  solution 
elusive.  I'herefore,  nun  h more  work  needs  to  he  done  not  only  with  the  Iheorv  hut 
al.so  with  experimental  data  gathering  hefore  a truly  valid  solution  will  exist  for  radar 
scattering  from  vegetation  type  areas. 


I.  Kiipirc,  I't  al.^  “Hcticrtion  arui  Tran<i|))ission  l)v  a Katidoin  Vlrdium,”  Hnilio  Srioncf^  Vol.  4.  No.  1 1.  Novcmttrr 
I'm.  pp.  lOfiT-iorr. 
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CONCl.USlONS 


'I'lif  r<>llt>\\ in;;  cotu'liisioii.s  ciiti  !«•  made  from  llic  study: 

a.  \ 111.  'ory  lias  bct'ii  d.'vt'lo|n‘d  lor  (’om|iuliii^  radar  liackscaltcr  Irom 
certain  ty|)es  id'  ve<:elalioii  by  usitii;  a v.-elor  renormalization  a|)|)roaeli. 

b.  No  riiforous  (|iiantilative  eoni|)arison  of  tlieorv  willi  e.\|ierinienl  was 
possible;  however,  (jiialilative  eom|iarisons  indicate  reasonable  awr.'emenl. 

e.  riie  theory  develo|)e.l  was  for  like  polarized  (llll  and  V V)  eomponenls 
only,  and  no  analytical  re.sults  were  obtained  for  the  depolarization  eomponenls. 

d.  \llboii.;li  no  explicit  solution  was  obtained  for  the  depolarization 
eomponents.  it  was  learned  that  one  eau.se  of  de|iolarizalioii  is  the  anisotropy  as.soeial.’d 
with  the  diele.  trie  fluetiialions  of  the  random  medium. 

e.  It  was  lound  that  the  mean  wave  in  the  vcfietation  attenuates  for  two 
different  reasons; 

I.  Average  lo.ss  due  l.j  water  eontenl  in  the  veijetation. 

I2.  Loss  due  to  multiple  .scattering. 

f.  A simple  quasi-qualitative  way  of  relating  the  statistical  parameters  of 
the  theory  to  the  actual  vegetation  parameters  was  provided. 

g.  Using  the  developed  theory,  it  was  found  that  the  baekseatter  coeffi- 
cient inerea.sed  the  fastest  tor  values  of  leaf  moisture  content  hetween  10  and  HO 
percent. 
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CAl.CULATiON  OF  TIIF,  FLFMFNTS  OK  TIIK  DiKLKCTRIC  TKXSOIl 

'I'hf  |>iir|><)Sf  (if  (his  :i|i|i(‘n(lix  is  (o  (Iclcrtiiine  (lie  lu-cdcd  clcincnls  of  llic  dicicc- 
trie  tensor  and  to  oldain  necessary  expressions  for  the  direction  cosines  of  the  effec- 
tive propagation  constant.  It  lias  already  been  shown  that  the  elements  and 
will  he  zero  for  an  i.sotrojiie  correlation  funetion.  The  elements  that  net'd  to  he  ( om 
pnted  are  c , e , e and  e^,.  We  will  consider  the  ( aleulation  of  e in  detail, 
and  only  state  the  results  for  the  other  elements. 

V ' 

1=  y'n(r-r')|f2(K)+  f3(R)}  e'- '*-■-’  ,ir' 

V ' 

We  let  B(r-r')  = I5(ll-£l)=  , (R  = Ir-r'l),  and  where  C is  the  correlation 

dislanec.  Also  we  will  change  the  variables  of  iiilcfjration  such  that  u = x-x'.  v = y-y’. 
and  w = z-z’. 

I=y J dwH(R)^f3(R)  + f3(R)}e-‘^-'-^ 

The  (geometry  of  the  K and  r - r'  vectors  is  given  lielow. 


w 


K 


Figure  I . Geometry  of  the  K and  r-r  Vectors 

will  now  rolalf  the  u.  w plane  until  the  direetion  of  the  new  w axis  (w')  eoineides 
with  the  direetion  of  K.  The  geometry  of  the  rotated  coordinate  system  is  given  helow. 


P 


riu'  intt‘i;ral  I now  Ix-coines 


= y^iu'y*  ds'j  aw'iHK)|f2(H)+ 


jK  -(r-r') 


If  we  eliange  now  to  the  spherical  coordinates  li.  6,  and  0,  the  intetcral  1 will  he: 

2n  n " 

1=  J' d(l)  J' do  J'  dK  K^sineiHK)  I f2(ll)  + (7^sin^0eos^0  + 

o o o 

^ iK.Kcos0 

2a7i‘t)s0sinOc()s0  + Of^ros^^)  t3(K)  / 

We  can  break  up  the  above  eoinplieated  integral  into  the  sum  of  a number  of  smaller 
integrals,  the  first  of  which  to  be  considered  is  Ij . 


1. 


dl{  sin01?(lf)  f2(R)  e' 


jK  R ros£> 


o o o 

15y  integrating  in  0,  we  have 


I,  = 


dlf  K^sinOlifK)  f2(K)c 

O 

Harrying  out  the  integration  in  0 will  yield 


jK  H ro.sC 


I , - dir 


/■"*  " 


kK 


R'  R(R)  f2(H) 


rh(’  stM’ond  inU’fjral  U>  l»t’  (■<)nsid(*nMl  llial  is  a part  of  I will  drsi^naU'd  I2  ant! 
{^ivcii  as  follows: 


2it  n 


jK  K cosO 


l2=J' d0 y* do  Jd[{  RHdn^O  B(R)7^  <)s^0f3  (R)e' 

0 O O 

W hen  the  integration  in  0 is  performed,  the  following  result  is  olitained  lor  I2  : 


73 


rr 

I2  = TT 

0 t) 

I'll.-  inU-nration  in  0 can  lu-  easily  |it'rl'omK‘<l  Uy  iisiii-;  (he  itiellii.d  of  iiilefiration  l.y 
parts.  W in’ll  lliis  is  coniplt’led.  I2  "ill  Im’i'oiiic: 


d|{l{^sill^0|{(l{)7^l3(l{) 


jK  U rosO 


4) 

riu-  third  integral  to  he  tdiisidercd  as  a part  of  1 will  he  di’si<;nalt’d  I3  and  is  •iim  ii 
lit'lovv  as 


I = / d0  / (10  / (IK  K^sin^  015(11)  2a7''<''^<^'»>'"/>f'3(K) ' 


jk  It  co.sO 


0 o o 


Intc'^raliin'  in  0 Ironi  zero  to  '2tt  will  <iiy(’  zero,  and  so  I3  is  identically  Z(  ro. 


13-0 


The  loiirth  intejiral  which  lortns  a part  of  I will  he  de.sifrnaled  I4  and  is  as  follows: 

jk  K <'<»0 


2 7T  IT 


1^  =J  ,10 J ,10 J (IK  K^sinOI5(K)a2cos20f3(K)e 

o o o 

(larryiiifi  out  the  inte<;ration  in  0 will  yield  the  (ollowiii};  result  lor  I4: 

rr  «" 

/•  JKKcosO 

14=  277  / (10  / (IK  K2sin015(K)a^„s^0f3{K)e 

()  O 

The  integration  in  0 can  he  carried  out  in  a straifilitforward  tiiainier  to  ftive 


(IK  15^15(15)  a^f3(K)  -2)  sin(KK)  + 2kKc„s  (kK)] 

4> 

.Since  I iscpial  to  the  sum  ,d'  I,.  I^.  I3  ^""1  '4,  write  an  (’(pialion  for  I inv(dNinc 

an  iiiiccration  in  K . 
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u 

^ l{^  -2)  sin(KK)  + 2Kllc<)s  (KK)]|iHK) 

licl'orf  proceedin';  lo  i:ile;;rale  in  |{.  we  f'liall  I'irfl  replace  K with  k,  wliieli  is  llie  first 

approximation  of  K wlien  = t).  It  eoidd  lie  ar;;iie(l  that  K should  hi'  kept  as  it  is  in 

the  integral  I,  and  then  the  solution  for  (and  all  the  other  elements  of  the  dielectric 

tensor)  will  he  a function  of  K.  This  would  make  the  final  solutions  for  K and  K 

. : ....  . yy  *'• 

implicit,  and  .some  type  of  iterative  solution  wiiuld  have  to  he  found.  However,  the 

ori;;inal  e(|uation  for  the  mean  wave  is  only  ^ood  lo  order  e^,  so  a eomplieated  iterative 

.solution  for  K and  K would  not  he  meaningful.  It  is  with  this  ar;ainienl  in  mind 

that  we  replace  K with  k in  the  inte;;ral. 

Jm  u>r,(ii)+  iai  r,(ii)  .,,«(kii)] 

(> 

+ [(k^H^-2)sin(kll)  + 2kHeo.s(kK)]J  IJ(U) 

When  the  approjiriate  (pianlilies  for  f2(R).  f3(li)  and  H(ll)  are  suhstituted  into  the 
ahove  e(|ualion  for  1.  and  after  some  al;;ehraie  manipulation,  the  followin';  re.sult  is 
ohiained  for  I: 


w hi' 


I 


drrk  ( I - )sin(kR)  + 47rj  - I ) 


<) 


sin(kR)  ^ 
R 


kR^  '[  kR 


+ 


i 

■i 


I .,,,.(1,11))]  . [).  kll,i,i(kli)  t ■„>s(kl(|)])  X 


k 


r 


liy  u.'iii>;  tlic  lollow in*;  iiilff;rals,  an  alp'lirait-  roiill  t an  lit-  tiltlainctl  lor  1. 

k 


sin  k\  il\  = 


k'  ^ a^ 


/' 

li 

(c) 

o 

/ ax  I /smloy  \ 

' r (“u — "''V 

O 


■'^"1  I \ I .a  I , 

— — -ctiskxl  (1\ - 1 - - t (.r‘  (a/k) 
kX  / k 


1 /sin  kx 
kx 


CDS  kx  j (lx  --  + a^/k^  ) cof’  (a/k)j 

() 

/ "•"  {-  ^ ■">  ^')}  '■'  = T If  * (; )’ 


icot 

-(t)  -(t)  -■ 

' (a/k) 

} 

^ r__ 

1 1 ^ i(:ia^  - 1) 

.4  - 1 / „ /L  \ 1. 

i-i(-' 

-7^) 

Lk^ 

+ a^J 

col  (^a/K)  ^ 

I’k^ 

('^ 

S)  (i)]  * 

(2a^  - 
k^ 

a 

' k 

cot-' 

;j7' 

2k^ 

[!*(!) 

'(j) 

■(1/ 

1 

voV 

■'  (a/k)] 

^ + (a/k)^  - i eof' 

(a/k) 

■Cif 

cor’ 

‘(I)] 

\n  a|(|ir(ixirnu(('  analytical  expression  for  e can  now  Ite  written  as  follow,- 


e,,  + k^  eM 

XX  O o s 
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Till-  iiiilv  (lilTinilly  in  ('oiii|iiiliii^  I ruiiii's  Iroin  lln-  air  ('tilani’ciil  li-rni  >ini'i‘  liolli  a ami 
k an-  i-i)in|)lt‘\.  I’lii-  liillow in;;  analysis  i'  |ir<i\ iili-ii  to  show  a solution  to  this  |irohli-ni: 


la-t  \=i-ot-‘  (a/k)  = (-ot-'  ( 


Lsiiifi  till-  tai  l that  k = - j q;^^.  wi-  havi- 


^ + +J^o)  <^o  +J“o) 

- l-ot r r 

(3  - ia 

J O J [_  O 


Writiii*;  tiu'  i‘(|nalion  in  tin-  I'orrii  olCotA.  we  will  have 

,.„1  \ = ^ = ( l/g  + Qq  ^ jPoXPo  ^ ja.))  - (i-j'^  -Hi-~-'’^)/2  _ j (f-^  + I ) 

siiiA  p^+a^,  (|J'  - I) 


j (•■-J + 1 ) (Pf^  + ^(l/e  + + jP^^ } (P^  + ja , ) (c-J-'  - I ) 

Solvin';  the  ahoye  i'(|nation  tor  in  trrins  of  8.  anil  P^^  will  f;iy(‘ 

\P^  -t-  21/3  I2(j3^  + a^)  + a (2p^  £ -i-  2a^  C + a ) I 

•>;a  1^0  J^o  '^o  o'  o ' o o' I 


P^  + 

~ M ri 


= r,  +jr2 


w lu'rr  r,  = I + 2 a C 

* O 


'^2=2P£ 


(ilianpn;;  from  rcclanpilar  to  |iolar  nooniinates.  \vi>  w ill  have 

,.2jA_j.  (J0o 


wfu-n-  r,,  = y r]+r\ 


0^  = tair'  (r2/ri)  vvlien  Pi  >0 

= TT  lair'  (r2/r, ) whi-n  r,  < 0 anil  r2  > h 
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0^^  = t;tir'  (ri/i'i ) - jr  wlirii  r|  < 0 ami  r2  < H 
riic  <|nanlilv  \ t an  now  he  ca.-ily  olilaimal, 

A = ‘,2  (0„  - jliii;,)  = . ol  ' (a/k) 

The  rest  ol'  llir  nonzero  eletnenl>  of  the  dieleelrie  (eiisor  can  lie  derivefl  in  llie 
same  manner  as  e was  derived.  InsU'ad  ol'  derivim'  eaeli  ol'  llie  oilier  elemenis  in 

XX  '■ 

delail.  we  shall  onlv  slale  llie  linal  residls. 


e ^ k^/k^  + k^  1. 

V y o a > 


where  1,  = 


iTT.P'  - i . ..r'  (»/k>. 


Up,,. 


or'(a/k) 


[j/d  Ma/k)^  - ^ eol-'  (a 


(il/k)  - col 


' (a/k)] 


-j^  ■ *]  ~ ^ j*' 


wlierc  I'  = 


t c/kl  + 

k + k' 


+ [a/k-(l  +aVk=)eol-'  (a/k)|  + 


M a, ] . |i 

’]  • i?  [5 


+ (a/k)^  - (a/k)  eof'  (a/k) 


- (a/k)^  eof'  (a/k) 


+ (a/k)^  - Y eof'  (a/k)  - (a/k)’  eof'  (a/k) 

K 


e k^  K 

X Z OS 


where  K = ^ eof-  (a/k)  - ^ 


k’ 


2k’ 


— + 


(a/k)^ 
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-(a/k)foC*  (a/k)  - (a/k)^  coC*  (a/k)j  + |i>/k 

-(1  + a"/k")c..|-'  (a/k)l  + [l  - ji  c.l-'  (a/k)] 

\Vf  will  now  coiisicler  an  approximate'  soliilion  lor  llu;  direction  cosines  a and  7. 
F(»r  a vertically  polarize'd  wave,  we  have 

K^^sini// = k^sinOj  where  sini//  = a; 

K = |3  - ja 

A XZ  '^XZ  J xz 

situ//  = (a  + jh)  sin0j=a 
where  5 = 

f)  = a k ) 

xz  o '^xz  xz' 

Wheiv  the  elements  of  the  dielectric  tensor  and  the  effective  pro|)agation  constant 
are-  being  ceemputed  anel  a value  for  a is  ne*e;eled,  we  can  simply  replace  hy 
its  first  approximatie)!)  whiedi  is  k.  'I'his  woulel  mean  that  anel  wejulel  he  re- 
plae-eel  hy  anel  re-s|)ee: lively  in  the  above  equations  feer  a anel  I*.  Alse)  when  7 is 
e-eernputeel,  we-  woulel  use  the  following  exfiression: 


7 = -ce)si//  = - y 1 - (a  + jh)^sin^0.  = - p,e  ■*  ' 
w here-  p^=  1 1 I - (a^  - h^  ) sin^  0;  P + da^  sin‘'  d.  \ ^ 


i tan-' 


2ahsin^Qi 


I - (a^  - 1>^)  sin^0j 


()ne-e  again  if  7 is  nee-eleel  to  compute  then  we  e-an  simply  replaea-  K hy  its  first 
appreexirnaliem,  k,  in  eereler  le)  obtain  a first  appreeximation  fe)r  7. 
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APPKNDIX  B 


EVALUATION  OF  AN  INTEGRAL 


The  integral  lo  he  evaluated  in  this  appendix  is  the  one  {jiven  hy  equation  (127), 
which  is  repeated  l)elow  : 


IT  2n 

^ L'  L'*  /*  I tt  /*  ] /*  I . ii-?  • -K/t*  -2jk  KMnxroh0sind 

< ^IM  > / «i0  H^'^inxe  t*  « 


O o o 


jqKcosx  ) * -jK(osx(k  +k.7)/2  ^ ♦ -jK  cosx(k' -k '*)/2 

" ^ +-\lS2‘'2‘‘  ^ 


* jKco.sx(k  -k;*)/2  ^ » 

+ A^2-Si‘^3^‘  +^2'^2‘=4‘'‘ 


jKcosx(k^+k  *)/2  ^ 


I he  integral  in  0 can  he  performed  directly  hy  using  the  following  result. 


/-jxcost 

c dt=2ir.l^(x) 


riie  expression  for<  E„|,  E*j,|  > then  becomes 


. (k  cos0.)T\ 

< F F > = ■ ” 


in  12  2 • jlKrosx 

(IK  K'‘smx  e 27r  e 


I /Ol  12  • ■ /]  \ i \ -jKco.sx(k'  + k'*)/2  * -jRcosx(k'  + 

,l„(2kj{smxsin0j)  j \^,  A^,,c,  e ^ ^ + SuS2^‘2‘' 


* jRcosx(k.^-k  )/2  * j 

^ S.\2'^3^‘  +A^2^\  c,C 


jRcosx(k^  + k^  )/2  I 


We  will  break  up  the  above  integral  into  four  separate  integrals,  the  first  id' which  we 
w ill  call  K , . 
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<*>  7T 

• r . r . „ jqK'osx  -jKco8x(k,'  + k'*)/2  , -K/ii 

K,  = 27r A^i  1 1 | / <1K  / dx  siiix  ,|^(2k^l{siiixsiii0j)  e e ^ ' ll^c 


O O 


let  u = cofix  siiix  = y I - cos^  x = }/  I - 

<iu  = - siiixdx 


ju|qK-|{(kAk  *)/2J  -li/e 
'■  R‘‘e 


K,  =27rA^,A‘,r,y’  dK  j du  .|„(2kJ{siM0, 

« 00-1  I 

K,  = 27rAy,  Ay,c,  ly  dl{  j du  \{^v  ]J2kJU\tiO.  /TiT)  cos  lu(qR-R(k^ -i- k'*)/2)] 

O -1 

«>  I 

+ j I dR  j du  R"  J„(2k^RsinO.  V I -u^  ) |ii  (qR-R  (k'  +k'l  )/2)| 
o -T 

The  second  term  on  the  ri<;litdiand  side  of  the  !d)ove  equation  is  zero  since  the  integral 
in  u is  from  -1  to  -t-l  over  an  odd  functioti. 


K,  =27rA^,A^,( 


I 

o -1 


du  R^  J^(2ky  Rsint?!^  V l-u^  I cos  (u(qR-R(k^-t-k  *)/2)l 


uw  I 

k,  =47rA^  , A*  ,c,  J dl{  J .lu  R^'-«/''  J^(2kJ{sin0i  /TiT" | cos  lu(qR-R(k'+k  J/2)] 

O O 

I'he  following  integral  can  he  used  to  evaluate  the  integral  in  u from  0 to  1. 

a 

/I  /I  2\  , . 1 _ •‘‘hi  (a  Vh^+  c^  J 

.1  (h  Va  - x^  ) cos  (cx)  (lx  - ' 

Vh^  + c^ 

() 

sin  [ V ^*^0  ^ I )/“l 


K,  - 47rA^  , A 


*,c,  y* (IK  K^  e-'‘/“ 


y 4k^ R^ sin^ Oj+R^  1 q-(k;-Hk'*)/2 ] ' 
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I,f(liii>;  l)^  = y ^ |(|  - (k^  + k'*  )/2| ^ 


K,  = 47rA 


(IK  U (■-“/‘^ 


siti  (I)  K) 


'I’lic  iiitfiiral  ill  K is  easily  e\aliiateil  liy  iisitii;  the  I'ollowiiif;  integral; 


sill  li\(lx 


2al. 


(a^  + 1)^ 


a>() 


K.  = 


II 


riie  second  inte<;ral  to  lie  considered  that  makes  n|)  < Kuu  K*||j  > will  he  desifinated 
K2  and  is  iiiven  helow  . 


K., 


“ n 


o o 


It  can  he  seen  that  K2  has  exactly  the  same  form  as  Kj,  and  therefore  the  result  for 
K2  <‘ari  he  written  down  directly. 


KttA  ,A*,,e„£^ 
1^  ^2  2 


|l+h? 


where  h,  = V 4k^sin^0j  + (q  - (kj  - k'*)/2]2 

The  third  integral  that  make.s  up  part  of  < E||,|  E*|||  > will  he  given  the  symhol  K3. 

«»  jx 

Kg  = 27rA*  j Ay2^‘3^^ J' dx  R^sinxe'^'''j^(2k^jKsinxsin(?j)  e"*^"^**  * e^  * * 

O O 

It  can  he  seen  that  K3  has  exactly  the  same  form  as  k, , and  the  result  for  kg  follows 
from  kj . 
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K..  = 


‘ 1 I + 1.^ 


wlicrc  1>2  = y + |<|  + (kj  - )/2|^ 


Tlif  liisl  iiitt‘>;r;il  that  inakf.s  ii|i  |turl  of  < F.jiu  K||(|  > will  he  dcsi^iiaU'il  K4  and  i.< 
irivfu  Ik'Iow  . 


00  TT 

K,j  = 27r.\^,,  A*^Cjy* ilU y* dx  K^di.x  JfsinxsinOj)  .■ 


jqKrosx  jK cosxtkj^ ■*  kj  )/J 


o o 


Afiain,  it  is  easily  seen  that  K4  has  the  satne  fortn  as  K,  and  so  the  result  can  simply  he 
written  as 


drr A A* e , / z 

I I l],2  e2  |^~  'vlH‘reh3=  Vdk^sin^O, + |q  + (k;+k;)/21 


'I'he  final  result  for<  K|||j  ^■(|((  ^ ean  now  lie  written. 


< K K > = 
^ '^'1111  *^'11  IK 


, .A  , A .,e,, 

I ^ yl  y2  i 


2k=«>\,.,WS,  ( A„a;,. 

^y2  ^ y I ‘■;t  , ^>2  ^y2*^'4  ) 

I + hj  8^ ) ^ j 


1 1 + h^  8"  P I I 


A 
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APPENDIX  C 


COMPUTER  PROGRAM  LISTING  FOR  THE 
CALCULATION  OF  THE  RADAR  BACKSCATTER  COEFFICIENTS 

P«0GRAH  RANDCHC INPUT, OUTPUT, TA PE5= INPUT, TAPE6=OUtPUTl 

RADAR  scattering  FRON  RANDOH  NEDIA 

USING  A VECTOR  RENORMALIZATION  FORMULATION 

COMPLEX  Qj,QK,UCaTI,QA,QLl,QLZ,UL3,QL4A,QL4,0LG,QL,QKl,QT,QKZP,Ql, 
IQZ, 03, 04, ns, 06,17,08, a ay  IN, OAY 10,0 ay  1,09, 010,01 1,012, 013, 014, 015, 
2aAT2N,QAY2D,aAYZ,QKZPC,aAYlC,aAYZC,QCl,QC2,aG3,QC4, 
3au0,Q31S,Q361,Q32S,C332,O33,OHHl,OHH2, QHH3,aHH4, QHM.QVVl, 
40VV2,QVV3,QvVh, OVV,QLLK,aALPH,  QGAM, 

5axxi,aXX2,axX31,QXX32, 0XX3,QXX41,0XX42,0XX4,QXX5,QXX6,QXX,QEXX, 
6QXZl,QXZ2,0XZ3,aXZ4, OX Z5 , QX Z ,0EX Z, OZ Z1 , QZ Z2 , OZZS , 0ZZ4 , QZZS , QZZ6 , 
ZOZZ,UFZZ,QXU,OXN,OX,OX2,OCaSC,OTN,uru,QHO, UNO, OFl, aF2,QF3,aF0,  026, 
SQ2Z,C2 8, 02 9, 01X1,030, 031, 032, 033, 034,01X2, 035, 036,037,038,039,040; 
904 1,042, 043, 0AZ1,0AZ2.0AY1,0AV2,0AV 1C, 0AV2C 
N=3 

DO  15  1=1, N 

REAOtG,  ICll  AER,SIG,Xi<aL,EPS 
131  FORHAT(4F10.Si 
pi=3.  lAisgzes" 

F=9.0E*9 

KRITEI6,25C)F, AER,SIG, XKOL.EPS 

250  FORMAT (IH  ,2KF  = , £11.5 , 5X , 4M AER=; FI 0 . 8, 5X ,4riSIG= , FIO 06, 5X ,3HKL= , FIB 
1.8,5X,4HEPS=,F10.8I 
E0=8.854E-12 
U=4.0EA7 
UO=PI*U 
H=2.*PI*F 

9a=W*SORT(UO*EO*AERI 

A1.Q  = SI&*30RI(UO»EO»A£R)/  (2.*E0»AtRI 

XLArt=3.*10*»(Sl 7F 

XKO=V(»SORTtUO»E01 

XL=XK0L/XK0 

OJ=(0. 0,1.0) 

C'K=BO-QJ*ALO 

RlN=BO**24ALO»(2.»XL*30*»2f2.*XL*ALO»*24ALO» 

Rl=RlN/iaO**2+ALO**2) 

R2=2.*30*XL 

R0=SuRItRl**2+R2**2l 

PHI0=ATAN2fR2,Rll 

QCOTI=0.5»(PHIO-QJ*ALOGIRO)I 

AlrREALIGCOTI) 

A2=AIMAG(QC0TI) 

0A=1./XL+QJ*0< 

A3=REALIQA) 

A4=AIMAG<0AI 

QLl=l./IQK**2+0A*»2r- 

A5=REALI0L1) 

A&=AIHAG  (OLD 

QL2=-OJ*OCOTIZaX»»T 

AZ=REAL(OLZ» 

Ad=AlHAGfQL2) 

OL3=-3.*aJ»IOAZOK-Il.+ tOArOKI**2)*OCOTI)/12.»OK*QKl 

A9=REALIQL3I 

Aia=AIMAG(QL3l 

QL4A=2./3.H0A/0KT*»2-IQA70K)*QCIITr” 

A11=REAL(0L4A) 

A12=AIMAG(0L4A> 

QL4=CQL4A-QCOTI*iaA/OX)»*3)/(2.’*QK*QlO 
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QL5=IQC0TI*Qft/a<-l.l/(QK*(JK| 

QL  = QL1»QL2<’QL3»QL4»»QL; 

QLLK=l.*QL»iXK0*XK0*EPS/aKI**2 

XB  = i<E  AL(QLLKI 

Yd=AIHAG(aLLK» 

PHI3=ATAN2IYB,XBI 
RH0a=SQRr<XB*»Z*Y8*»2Y 
QKl=QK*SQftHRHOB)  •CEXPiaj'PHIB/Z.I 
B1  = R£AI.IQK1) 

AL1=-AIMAG(QK1I 

AA=XK0*B1/IB1**Z*AL1»*Z> 

eB=XK0*ALl/{Bl**Z*ALl»*2l 

THETA=0.0 

HRirE<6tllU 

111  FORMAT IZX> BHTHETA, lAX, IHP.lOX, ZHPl f SX ,?HSIGHH, 11 X t 5HSI GV V. 10 X, 3HXK 
lI,lQX,BHSIGCMt 
lOZ  T=PI»THETA/180. 

SM=SlN(ri 

CN=COS(TI 

Y1=Z.*AA*B3»SN*SN 

X1  = 1.HB3**Z-AA*»Z»*SN*SN 

RH01=lXl»*2*Yl»*ZI»*0.Zf 

PHI1=0.5*ATAN2(Y1,XH 

P=RH01*<Bl*SIK(PHri) ♦AL1»C0S«PHI11I 

Q=RH01*I31*COS(PHI1)-AL1*SIN{PHI1J7 

QT  = 2.*XK0*CN/ia-QJ»P«'XK0*CN) 

XKX=XKO*SN 

XKY=0.0 

XKZ=SgRTtXK0*»2-XKX**2-XKY»»2» 

YZ=-2.*ALO*BO 

X2=80**Z-ALO*»Z-XKX*»2-XKY**Z 

RriOZ=SQRT<Y2**Z*X2**2» 

PHIZ=ATAN2(Y2,X2I 

!JKZP=SQRriRHOZ»*CEXP«aj*PHIZ/Zi» 

QKZPC=CONJG(a<ZP) 

XKR=SQRT tRHGZ)*C0SCPHIZ/2.J 
XKI=SaRTlRHSZJ»SINtPHIZ/2.T- 

Ql=QKZP*X<Z*»2*a<ZP*XKX»»2+XKZ*QKZP»»Z«-XKZ*XKX**2 
Q2  = XKZ»IXKZ*Q<ZP) *EP3*3T*tXK0*XKX»XKY»  **2/ (Z.’QJ'OKZP) 
Q3=ai*XKZ*QKZP*XK0*XK0*EPS»gr/ (2.*QJI 

Qfc=XKZ*ai*f X<Y»*2+QKZP**2) • X KO *XKO •EPS*QT/ (2. • QJ*QKZPI 

Q5=-GKZP»IQKZP»XKZJ»(XKX»XKY>**2 

Q6=XKZ*CXKY**2+QKZP*»Z1»Q1 

Q7=-XKZ*<XKZ»QKZPJ* (XKX*XKYI**2 

Q8  = gi*IQKZP*X.<Y**2*QKZP*XXZ**2l 

QAYlN=!32*a3-Q«» 

QAYlC=Q5  + Q6*-QZ  + a8 
QAYl=aAYlN/QAYlD 

Q9=XKZ»Ql*fX:<T**2*QKZP»»2)*XK0»XIC3»EPS'»aT/lZ.»QJ»QKZPr 

Q10=-X<Z»IQ<ZP*X<ZJ ’iPS'ClT* {XKX*X<Y*XK01**Z/ (2.»aJ*QKZP) 

ail=Ql»XICZ*UKZP»XKO*XKO*EPS*ar/(Z.*QJl 

Q12=-QKZP*<0KZP*XKZ» ' (XKX’XKYI **2 

Q13=XKZ»  «XKY**2K1KZP**2)  *111 

Q14=-XKZ*IXKZ»aKZPJ*(XKX*XKy|**2 

ai5  = QKZP'*IXKY*»2tXKZ»*2)*ai 

QAYZN=Q9»Q10»ail 

QAY20=Q12*Q13«'Q14»Q15 

QAY2=QAY2N/QAr20 
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QAYlC=CONjC(aAiril 
QAy2C=C0NJG(aAr2» 
0C1=1./(?.»P-QJ»QKZP*UJ*QKZPC) 
QC2=1./(2.»P-QJ*QKZP-4J*QKZPCI 
QC3=l./(2.*F  + aj»aKZP«-QJ*QKZPC) 

QC4=i./  (2,»p*aj»aK.zp-aj*UKZPci 
Q30= t2.»XK0»SU) »*2*{Q- (QKZP*QKZPCI/2.)*»2 
8Bas=REALIQB0l 
bao=saRT (oBOsi 

Q31S=(2.*XK0*SNI**2*tQ-(QKZP-QKZPCI/2.l*»2 
XlB=PEAHaBlS) 
yia=AIMAG<CBlS) 

RHOai  = SQRT  {X13*»2+T13”2) 
PHI31=ATAN2iyi3,XlB» 

QBBl=SQRrtRHaBl) ♦CEXP<QJ*PHIB1/2.1 
Q32S=(Z . ♦XK0*SN) **2 ♦ t a+ ( QKZP-QKZPC) /Z. ) 
X2B=REAL(Q32S) 
y23=AIHAG<Qa2SI 
P.HQa2  = SaRT(X2B**2»y  23**21 
PHI32=ATAN2 (y20, X2BJ 
a332=SQRTCRH0a2)*CEX3 (aJ*PHIB2/2.l 
aB3=l2.*XK0*S(4I  •*2+(3*  taKZP*aKZPCI  '2.1 
33=REAL(QB3I 
3a3=SQRT(B3l 

QHHl  = QAyi»QAYir:*aClZ(i.*(aB0*XU1'»*2)**2 
GHH2=-ljATl*(iAy2G»aC2/  J 1 . + ! SBBl*  XL  I **2  .*  **2 
QHH3  = uAy  LC»QAy2»QC3/  (!.♦  ( C1B82* XL  » * *2» **2 
UHHi*  = QAy2*QAy2C*aC‘*Z<l.  + CBB3*XL)**2l**Z 
t;HH  = QHHl*aHH2PaHH3+QHH4 
HH=REAL(QHH) 

SHH=8.*IXL»*3)*IXK0»:N»**2»IHH1 
QaS4H=lG .'ALOGIQ (SHHI 
QALPH  = XKO*SN/(30-(1J*ALOI 
AHl=B0*XKC/(33**2+AL0**21 
BHl=AL0*XK0/ia0**2*AL0**21 
yil=2.*AHl*BHl*SN**2 
X11=1.-IAH1**2-3H1**2) *SN**2 
RH021=(Xll**2*yil**2)**0.25 
PH121  = 0.5*ATA‘J2<yil,Xlll 
QGAM=-RH0  21*CEXPI-QJ»'>HI21» 
QXXl=aGAM**2/(GK**2+QA**2l 
QXX2=GJ*(3.*QALPH**2-l.)*aC0TI/tQK*,GKI 
QXX31=QJ*3.*t2.*QnLP3»»2-QGAM**2J/ (2.*QK*TJK1 
QXX32  = uA/QK-(l.  + (GAZat<»**2)  *00011 
QXX3=aXX31*QXX32 

QXX41=I2.*QALPH**2-QGAH**2I7<QK*QKI 
QXX42=l.-(QA/a<)*QC0II 

QXX4=QXX41*GXX42 
GXX5=aL4*3.*QGAN**2 
QXX6=-2.*QL4 

QXX=QXXl*axX2*3XX3f0XX4*QXX5*aXX6 
QEXX=(QK/XKO»  •* 2*axx»  J XK  0*EPS) 
QXZl=3.*aJ*QGAM*aALPH*aC0Il/CQK*QKI 
I QXZ2=-3.*QGAM*aALPH*aL4 

; GXZ3=-QGAM*QALPH*QL1 

QXZ4  = 9.*0J*aGA>1*aALPH*aXX32/(2.*QK*QKI 
, 0XZ5=3.*QGAH*QALPH*UXX42/(aK*aKI 

I UXZ=QXZl*QXZ2*aXZ3*QXZ4*QXZ5 


Q£XZ=QX2»IXK0»EI»Sl»*: 

QZZ1=QL1»(QALPH)»*2 

QZZ2=aj*l3.*9GAi**Z-l. ) *00011/ «QK*QK1 

QZZ3  = QJ*3.*«2.*aGAH**2-QALPH**Zl  •QXX32/ < 2. ’aK'dlCT 

UZZ<»=  (2.*a&AH**2-aALPH*»2)  'OXXifZ/ ( QK*QK) 

ClZZ5=3.*0ALPH»(JALPH*aL4 

QZZ6=-2.»QL‘* 

QZZ=QZZl»QZZ2»aiZ3»QZZ4+QZZ5*(lZZ& 

UEZZ=  (QK/XKOI  •‘■Z  + QZZ*  tXKO’EPSJ  **2 

QXD=UEXX*QALPH**2*QEZZ*QGAM**2*2.*aGAM*QAcPH*QEXZ 

QXN=aEXX*QEZZ-UEXZ*»2 

QX=QXN/QXO 

X=REAL<QX1 

r=AIMAG»(3X» 

RH0  = SaRHX**2*r**2l 
FHI=ArAN2(rtX) 

UK2=XK0*SQkT(RH0l*CEXPCQJ*PHI/2.» 

aALPH=XK0*SN/’(lX2 

BETA2=REAHQK2) 

ALPH2=-AIMA&(QK2) 

AH=SETA2*XK0/(aETA2**2+ALPH2**2l 

8H= ALPH2*XK0/(aETA2*»2+ALPH2**2l 

RH0  2X=l.-<AH»*2-dH**2»  *SN**2 

RH02Y=2.* AH*3H»SN*SN 

RH02=CRHaZX»»2+RH0ZY**2) ••0.25 

PHIZ=0.5»ATAN2 (RH02Y, RH02XI 

QC0SC=RH02*CtXP«-QJ*PHI2l 

QGAM=-QCOSC 

2TS’  = 2.*aK2*XK0»CM 

QrD=QK2*XKD*CN+aC0SC»XK0*»2 

QT=aTN/QrO 

P1=RH02» (ALPH2»OOS(PHI2)*BETA2*SIN<PHI2)) 

Q=RH02* {BETA2»C0S(PHI2i -ALPK2* SIN t?HI2) } 
QM0=EPS»M»U0*QT*RH02»CEXP«-aJ*PH121 •XKO*XKO/ai02 
QN0=XK0*»3»EPS»R*Ua*aT*SN/ (QK2**2) 

QF1=<QKZP*XKZ> •XKY*»2*QKZP*XKZ**2»XKZ*QKZP*»2 
QF2=XKX*XKY*taKZP»XKZ) 

QF3=XKZ*(XKZ*QKZP»QKZP»*2l*IXKZ*aKZP)*XKX**2 

QF0=QF1*QF3-QF2**2 

Q26=XKX*XXZ*aN0»aFl/(2.*Q3> 

Q2Z=XKX**2*XKZ*aMa»0Fl/(Z.»aJ*aKZP) 

Q28=-XKY*XKZ*aN0*QF2/C2.*QJI 

Q29=-XKY*XKX*XKZ*QI10*aF27I2.»aJ’QIC2P3 

QAXl=IQ26*Q2Zf a28»Q29»/QF0 

Q30=-XKX*XKZ*aNO*QFl/ (2,*aJI 

a31=QJ»QM0*XKZ*QKZP*QFl 

Q32=-XKX**2»XKZ»Q10*aFl/ (2.*QJ*QKZPI 

a33=XKY*XK2*QN0*aF2/(2.*QJI 

Q3i»=XKY»XKX»XKZ»QH0*aF2/(2.*QJ»QKZP) 

OAXZ=«a30»Q31*a32  + Q33  + Q34l/'QFO 

Q35=XKY*XKZ»aN0»aF3/(2.*QJ» 

Q36=XKY*XKX*XKZ*aM0*aF3712.*QJ*QKZPI 

a37=-XKX*XKZ*QN0»QFZ/l2.*aJI 

038=-XKX**Z*XKZ*QMC*QF2/<2.*QJ*QKZP» 

QAYl=IQ35*Q36»a37*Q38) /QFD 

C39=-XKY*XKZ»QN0»QF3/ (2.»aJ) 

Q4  0 = -XKY*XKX»XKZ*Q'10»aF3/(2.*aJ*nx7P» 
a%i=xKx*xKZ*aNo*aFr/f 2.*QJi 


87 


Qi*2  = -QJ»Qno»XKZ*QKZP»QF2 

Qi.3  = XKX*»2*XKZ»QM0»CFZ/i2.»QJ»UKZPI 

QAY2=  (Q39f  Q<«0»Q‘>l»Qi>2*Q43)/QFO 

Q4Zl  = XKX*aAXl/XKZ*XKY»ClAYlZXKZ 

aAZ2=XKX*aAX2/X<2»XKY*QAY2ZXKZ 

UAYl=-<aAxl*CN*aAZl*SNI 

aAV2=-(QAX2*CN*3AZ2»SM» 

QAVlC=CONJGCaAVll 

GAi/2C=C0NJG((lAI/2» 

aci=i./<2.*pi-aj*aKZP»Qj*QKZPci 

a:2=i./(2.*pi-aj»aKZ»-Qj»a<zpc» 


QC3=1./<2.*P1»GJ 
004=1. /(2.»P1+QJ 
QBJ=(2.»XK0*SN) 
BBQS=REALfOBOI 
BaO=SQRT (BBOSl 
QB1S= I2.»XK0»SNI 
X13  = P,EAL(Q81SI 
YlB=AIMAGtaBlS’ 
RHObl=SQRT{Xia*» 
“HI J1=ATAN2IY13, 
QBB1=SQRTIRH0B1> 
Qa2S= I2.»XK0*SH) 
X2fl=REAHQB2S> 
yZB=AIMAGIQB2S) 
“HJ32=S0Rr<X2e*» 
PHI32=AIAN2 (y23, 
aaa2=saRTtRH032) 
a33=;2.*XX3»3N) * 


*aKZ“»oj*a<zpci 

•QKZ“-aj*a<zpci 

»2»<a-(QiCZPtQKZPC»/2.l»»2 


**2*(a-<QKZP-QKZPC)/2. »**2 


2*Yia**2l 

XIB) 

•CEXPIQJ»PHIB1/2.1 

»*2*ia+«aKZP-aKZPc>/2.)»*2 


2»Y23»»2I 

X2B} 

*CEX“iaj»PHIB2F2i1 
* c* (0* I QKZP  V QK  ZPCi/Z.)'*2 


B3=REALia33) 

Ba3=soPriB3) 

OVYl  = aAVl»aAtflC»aCl/<i .♦  I330*XL***2I**2 
QVtf2=QAYl*QAY2C»aC2/ H.*IQB81*XL)*»2»**2 
G\/V3  = 0AVlC*aAV2*aC3/«l.  + IQ3B2»XLI**2»»»Z 
QVY4=QAV2*aAW2C*aC4/  II  .HB33*XL)  **21  **2 
Q7W  = 0\/Vl  + 0>/V2*-atfV3  + QVV4 
V7=P.EAL{QVV) 

SVV=8.*(XL*»3l*<H»E0*CN>**2*Hf») 

OBSVV=1I).*ALOG10  (SVVI 
NK=THET A 
IF(NKI15.260»2&1 

260  SIGCO=SHH 
OBSCM=10.*AUOG10ISIGCtn 
GO  TO  16 

261  SIGCC=SIGCO*CN 
OaSCH=10 .'ALOGID ISIGCCI 

16  WRITE (b.lOeiTHETA.P.Pl.DaSHH.OaSVX.XKIfOBSCI 

106  FORMAT CIH  » 2X,F 5. 1 , lOX , F 7. 3» 5X , F Z. 3, 3X, F? . 3, 11 X, FT .3 ,5X, FS. 3 ,5X, FT 
1.3» 

14  IFITHETA-BO.I 107,15,15 

107  THETA=THETA+10. 

GO  TO  102 

15  CONTINUE 
STOP 
ENO- 
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rclu(i\c  (iifli-ctrir  coiistuiil  of  tlii'  raiiiloiii  ni«‘iliiiiii 
avtTuuc  relative  ilieleetrie  eon^laiit  of  llie  raiitloiii  iiieiliiiiii 
average  eoinlnelivilv  of  the  raihloiii  ineiliiiin 
slaiitlard  ilevialioii  of  llie  dieleetrie  fliieliialioiiN 

random  portion  of  tlie  relative  dieleetrie  eon.-tant  in  tlie  random 
medium 

autoeorrelation  funetion  of  fiir) 

ra<lian  fretpienev  of  tlie  ineident  wave  | 

relative  mean  eomplex  ilieleetrie  eonstant  of  a fore>t  of  leave^  ^ 

volume  of  one  leaf  ! 

total  numlier  of  leaven  in  tlie  forent  volume  \ 

relative  eomplex  dieleetrie  eonstant  of  one  leaf 
volume  of  the  forest  wliieli  is  solelv  air 
dieleetrie  eonstant  of  air 
total  volume  of  the  forest 
(lirlrctnc  roiislaiit  of  free  spaee 
anffle  of  ineidenee 

total  eleetrie  field  in  the  random  medium 
total  mut;netie  field  in  the  random  medium 
del  operator 

permeahility  of  free  spaee 
scattered  eleetrie  field  in  tin-  random  medium 
infinite  spaee  dyadic  (Ireen’s  function 
propajiation  constant  in  free  spaee 

mean  hack^roiind  propagation  constant  in  the  random  medium 
unit  dyadic 

mean  eleetrie  field  in  the  random  medium 

effective  propagation  vector  of  the  mean  wavi-  in  the  random 
medium 

relative  complex  dielectric  tensor 
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(‘I<  iiu‘iils  of  llic  r<-lali\i‘  ilirli't  lrit'  IriiMtr 

unit  vrriors  in  tin-  > aiut  /. 
ilim'liiin  i-itsiiio  ttf  tli<‘  K vrrinr 
rfflcrlion  ctM'fririnil  (liori/.onlul  |i«ilariy.ali>Hi) 
traiiMiiissioii  rtN'fficH'iil  (Intrixonlal  ixdari/.alioii) 

atiriiiialion  of  (In-  iiH-an  ua\i‘  in  llic  rainloin  nii'ilinni  (Intri/.oii- 
(al  |i«(lariyalii>n) 

Konrirr  \arialil«‘.s 

srallrn‘<l  <‘l<‘«'lri«-  ndil  in  Ihr  u|t|MT  inrdinni  (air) 

radar  ltu«-k.>a-allfr  ctM'frii'H'iil  for  lli«'  caM*  of  liori/.onlal  |tolari/a- 
lion  Iran.Miiil  — liori/.ontal  polari/.alion  r<‘t-i-i\<* 

rcfli'i'lion  ( (M-fricicnl  (vcriiral  polarir.aliiMi) 

Iransniission  riH-fniU'iil  (vertical  polarixalion ) 

altcniialion  of  llic  inean  wave  in  llic  randinn  niedinni  (verlii  al 
p(dari/.alion) 

radar  liaek.M  ulter  co«-ffi<  i«‘nl  for  the  easa*  of  vertical  |iolari/.ati<Mi 
Iransinil  — vertical  polari/.atiun  n'ceive 
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